Master crossover behavior of parachor correlations for one-component fluids. 
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The master asymptotic behavior of the usual parachor correlations, expressing surface tension a 
as a power law of the density difference pl — pv between coexisting liquid and vapor, is analyzed 
for a series of pure compounds close to their liquid-vapor critical point, using only four critical 
parameters (/3 C ) _1 , ot c , Z c and Y c , for each fluid. This is accomplished by the scale dilatation 
method of the fluid variables where, in addition to the energy unit (/3 C ) _1 and the length unit 
a c , the dimensionless numbers Z c and Y c are the characteristic scale factors of the ordering field 
along the critical isotherm, and the temperature field along the critical isochor, respectively. The 
scale dilatation method is then formally analogous to the basic system-dependent formulation of the 
renormalization theory. Accounting for the hyperscaling law fcl = we show that the Ising-like 
asymptotic value ir a of the parachor exponent is unequivocally linked to the critical exponents n, 
or S, by -j^j — d _ l ^_ v - ) — (here d — 3 is the space dimension). Such mixed hyperscaling laws 
combine either the exponent n, or the exponent 5 which characterizes bulk critical properties of d 
dimension along the critical isotherm or exactly at the critical point, with the parachor exponent 
ir a which characterizes interfacial properties of d — 1 dimension in the non-homogeneous domain. 

Then we show that the asymptotic (symmetric) power law (a c ) d_1 [3 c a = Dp (^ PL ^ p Pv j is the 

two-dimensional critical equation of state of the liquid-gas interface between the two-phase system 
at constant total (critical) density p c . This power law complements the asymptotic (antisymmetric) 
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of the three-dimensional critical equation of state for a fluid of 

density p 7^ p c and pressure p 7^ p c , maintained at constant (critical) temperature T = T c [p p (p p , c ) 
is the specific (critical) chemical potential; p c is the critical pressure; T c is the critical temperature]. 

d 

We demonstrate the existence of the related universal amplitude combination Dp [Dp J = Rdo- = 
universal constant, constructed with the amplitudes Dp and Dp, separating then the respective 
contributions of each scale factor Y c and Z c , characteristic of each thermodynamic path, i.e., the 
critical isochore and the critical isotherm (or the critical point), respectively. The main consequences 
of these theoretical estimations are discussed in the light of engineering applications and process 
simulations where parachor correlations constitute one of the most practical method for estimating 
surface tension from density and capillary rise measurements. 

PACS numbers: 64.60.Ak.;05.10.Cc.;05.70. Jk;65.20.+w 



INTRODUCTION 



Most of the phenomenological approaches for mod- 
elling the fluid properties in engineering applications are 
commonly based on the extended corresponding-states 
principle (author?) [1]. In this scheme, the estimation 
of thermodynamic properties can be made using multi- 
parameter equations of state, that account for increasing 
molecular complexity by increasing the number of ad- 
justable parameters. Such engineering equations of state 
(whose mathematical forms must be compatible for prac- 
tical use in fluid mixture cases), are then generally con- 
venient tools to estimate a single phase property with 



sufficient accuracy (author?) [1]. However the knowl- 
edge of properties in the nonhomogeneous domain (au- 
thor?) [2], such as the surface tension a, the capillary 
length lea. the density difference Ap^v = Pl — Pv be- 
tween the coexisting liquid and vapor phases of respec- 
tive density and pv , are also of prime importance 
to gain confidence in fluid modelling and process sim- 
ulations (geological fluid flows, assisted recovery of oil, 
storage of green house gases, pool boiling phenomena, 
microfluidic devices based on wetting phenomena, etc). 
Therefore, a large number of related phenomenological 
laws, referred to as ancillary equations, have been pro- 
posed in the literature (author?) [1, 2, 3] to calculate 
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such properties in the nonhomogencous domain. This 
complementary approach generally leads to unsolvable 
mathematical differences with values calculated from the 
equations of state, and increases in a substantial man- 
ner the number of adjustable parameters to account for 
complex molecular fluids. 

Before focussing on the specific form of ancillary equa- 
tions between a and Ap LV (author?) [4], the so-called 
parachor correlations (author?) [1, 5], it is interesting to 
recall the two well-known practical interests of the fluid 
modelling based on the extended corresponding-states 
principle: i) the thermodynamics properties of a selected 
pure fluid are fully specified from a few fluid-dependent 
parameters such as for example its critical coordinates T c 
(critical temperature), p c (critical pressure) and v c (crit- 
ical molar volume) in the original and simplest form of 
the corresponding-states principle (author?) [6]; ii) the 
most convenient tools to estimate the fluid phase sur- 
face, including the two-phase equilibrium lines, are pro- 
vided by the cubic and generalized van der Waals equa- 
tions (author?) [1]. Our main objective in this intro- 
ductive discussion is then only to recall the number and 
the nature of the most usefull macroscopic parameters 
used in such engeneering equations of state (for a review 
see Ref. (author?) [7]). For more detailed presentations 
of the basic understanding from a rigorous microscopic 
approach of the molecular interaction and the theoreti- 
cal background for devclopping better functional forms 
of the pressure-volume relationship see for example the 
Refs. (author?) [1, 8, 9, 10] and the review of Ref. 
(author?) [11]. 

It was well-known (author?) [6] that only the in- 
ert gases (Ar, Kr, Xe) can obey the two-parameter 
corresponding-states principle (i.e., an energy unit and 
a length unit mandatorily needed to compare dimension- 
less thermodynamic states for same values of the dimen- 
sionless independent variables, admitting that the molar 
mass of each one-component fluid is known). This re- 
strictive conclusion was founded on results obtained in 
building unique functions of the reduced thermodynamic 
variables, examinating many thermodynamics proper- 
ties, such as the density difference between coexisting 
liquid and vapour phases, the saturated vapour pressure 
curve, the second virial coefficient, etc. 

Considering the modelling based on statistical me- 
chanics (author?) [8, 9, 10], this two-parameter 
corresponding-states description can be validated from 
the restrictive compounds made of spherical atoms 
with centro-symmetrical forces (such as precisely the 
inert gases mentioned above). The short-ranged 
space (r) dependence of intermolecular pair potentials 
u (r) (author?) [12] can be written for example in the 

form u (r) = eL.i^ LJ (577) wriCre ^ LJ i s the Lcnnard 
Jones (12-6) universal function (author?) [8]. The two 
quantities clj and ctlj are scaling (energy and length) 
parameters which characterize a particular substance. 
Compounds which obey this kind of universal potential 



function with two microscopic scaling parameters are said 
to be conformal (author?) [1, 9, 10]. 

On the other hand, some intermolecular potential mod- 
els with attractive interaction forces of infinite range, 
have given physical reality to the famous form (cubic 
with respect to volume) of the van der Waals (vdW) 
equation of state (author?) [13], separating then the 
repulsive and attractive contribution to the pressure- 
volume relationship estimated from the generalized van 
der Waals theory (author?) [7]. Although the two pres- 
sure terms of the original van der Waals equation do 
not quantitatively represent the true repulsive and at- 
tractive forces, the introduction of the two characteris- 
tic constants for each fluid - its actual covolume 6, not 
available to molecular motion due to a finite diameter of 
each repulsive molecule, and the amplitude a of the pres- 
sure decrease due to the intermolecular attraction -, has 
proven to be extremely valuable for the representation of 
its properties. Thus, after expressing the values of the 
van der Waals parameters a and b at the critical point, 
the unique function 2- = f vdw |l of the original 

van der Waals equation conforms to the two-parameter 
corresponding-states principle since v c depends unequiv- 
ocally on T c and p c , through the unique value of the 
critical compression factor Z^ w = = |. 

So that, at the macroscopic level, practical formula- 
tions of the two-parameter corresponding-states princi- 
ple employ as scaling parameters the critical tempera- 
ture T c (providing energy unit by introducing the Boltz- 
man factor fcs), and the critical pressure p c (providing a 

length unit through the quantity ( fc ^ Te ^ expressed for 

space dimension d = 3), and seek to represent thermo- 
dynamic properties, thermodynamic potentials and re- 
lated equations of state as universal (i.e., unique) di- 
mensionless functions of the new reduced variables 5- 
and (or ^-). However, although this principle only 
applies to conformal fluids, it is easy to show that it 
allways generates unreductible difficulties to obtain sat- 
isfactorily agreement between theoretical modelling and 
experimental results, especially for the two-phase surface 
approaching the liquid-gas critical point. For example, 
the potential parameters elj and ctlj of a Lennard- Jones 
(12-6) fluid evaluated from different thermodynamic and 
transport properties of the same real fluid tend to be sig- 
nificantly different than the ones directly obtained from 
their relations to the critical point coordinates (although, 
according to the molecular theory, the calculated critical 
compression factor Z\l ~ 0.290 remains the same for 
all these conformal fluids). Moreover, real atoms like 
Ar, Kr, and Xc are definitively not conformal (Z c , a c .u, 
and then Y c , etc., are not strictly constant numbers 
(author?) [9]). Similarly, the well-known breacking up 
of van der Waals equation of state occurs immediately, 
noting that the value Zq£ w = 0.375 significantly differs 
from the Z c values of real fluids [ranging for example 
from Z c (H 2 0) ~ 0.22 to Z c ( 4 Hc) ~ 0.30], especially in 
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the inert gase case [for example Z c (Xc) = 0.286]. 

As stated above in the developpement of its simplest 
form from fluid state theories, a two-parameter descrip- 
tion does not hold for real atoms and a fortiori for 
molecules with more complex shapes and interactions. 
Indeed, for compounds with non associating and non (or 
weakly) polar interactions of nonsphcrical molecules, also 
referred to as normal compounds, the deviations from 
the two-parameter corresponding-states modelling were 
most often described by one additional parameter, the 

so-called acentric factor, uo = — 1 — log Psatl ~ — — , 

Pc 

proposed by Pitzer (author?) [14]. The acentric fac- 
tor was defined from the reduced value p* = of the 

Pc 

saturated vapor pressure p sa ± (T) at the reduced value 
T* = ^- — 0.7 of the vapor saturation temperature, such 
that it is essentially ui ~ for inert gases Ar, Kr, and Xc. 
An inert vapor condensates at one tenth of the critical 
pressure at T* = 0.7, while a vapour of more complex 
molecules condensates at lower relative pressure, leading 
to cu positive, so that, the larger and more elongated the 
molecule, the larger to, due to an increasing contribution 
of the attractive molecular interaction. Thermodynam- 
ics properties of the normal compounds can then be de- 
scribed by unique functions of the three parameters T c , 
Pc w. 

It is also well-established (author?) [1, 7] that this 
three-parameter coresponding state modelling can be ac- 
counted for by using a three-parameter equation of state, 
for example the frequently referred Soave-Rcdlich-Kwong 
(author?) [15, 16], Peng-Robinson (author?) [17], or 
Patel-Teja (author?) [18, 19] cubic equations of state. 
When a third parameter is introduced into a cubic equa- 
tion of state, the critical compression factor Zq c becomes 
fluid-dependent, as stated for real fluids. Unfortunately, 
although a three-parameter equation can be forced to the 
correct Z c , only better overall improvement of the accu- 
racy in estimations of the phase surface is obtained when 
its calculated value is greater than the real one (for ex- 
ample Zf c K = 0.333 and Z% C R = 0.307 for Redlich-Kwong 
(RK) and Peng-Robinson (PR) equations of state, re- 
spectively, while Patel-Teja (PT) equation of state treats 
this calculated compression factor Zq c t as an adjustable 
parameter). More generally, such quantitatively inaccu- 
rate calculations are due to the relative rigidity of the cu- 
bic form (which limits the quality of the representation of 
derivative properties), added to fundamental limitations 
of analytic equations close to the critical point (which 
generate mean field behaviors of fluid properties). 

Moreover, useful precise measurement of the satura- 
tion pressure curve p sa t (T), introduces a critical lim- 



iting (dimcnsionlcss) slope ol Ci r 



aiog(T) 



T=T C 



9p s 



dT 



T=T C 



at T c , as another fluid characteristic pa- 



rameter, also known as the Riedel factor (author?) [20]. 
Anticipating then result of the next Section which in- 
troduces the critical number Y c = y- 



T=T C 



[see Eq. (19)], we note the relation a c ,R = Y c + 1 be- 
tween the Riedel factor and the dimcnsionlcss number 
Y c . Therefore, as an immediate consequence of the real 
location of the liquid-gas critical point in the experimen- 
tal p,v,T phase surface, the addition of Z c and a c ,R, 
(or Y c ) to uj appears as a useful parameter set incre- 
ment, able to describe deviations from the two-parameter 
corresponding-states principle based on T c , p c knowl- 
edge. Obviously, any three-parameter corresponding- 
states modelling needs implicit dependence between to, 
Z c , and u c .r (or Y c ), which provides bases for a large 
number of three-parameter corresponding-states models 
by developping empirical combinations between u>, Z c , 
and a c ,R (or Y c ), such as ones where Z c is linearly re- 
lated to cj (author?) [21] for normal compounds. This 
constrained situation to reproduce the critical point loca- 
tion in the p, T diagram, is certainly the most important 
practical reason why the three-parameter cubic equations 
of state (which allow a fair thermodynamic description 
of the normal compounds (author?) [1], including their 
intcrfacial properties (author?) [22]), are the most pop- 
ular equations of state developped again today for indus- 
trial process design. 

However, the three-parameter corresponding-states 
modelling still remains not appropriate for describing 
highly polar and "associating" fluids (such as water or 
alcohols for example). At least an additional fourth pa- 
rameter is needed to extend the corresponding-states ap- 
proaches, which leads to multiple routes to account for 
this increasing complexity of the microscopic molecu- 
lar interaction. Several empirical expressions have been 
proposed for this fourth parameter increment, such as 
the one introducing the Stiel polar factor (author?) 
[23] for example. Again a myriad of four-parameter 
corresponding-states models can then be defined using 
T c , p c , and practical combinations which provide only 
two independent dimcnsionlcss numbers chosen among 
the critical compression factor, the Pitzer acentric factor, 
the Riedel factor, the Stiel polar factor, etc. For example, 
Xiang (author?) [3], noticing that polar and nonpolar 
coumpounds may have similar u>, but different Z Cl or, in 
other words, that the relation between Z c and w, that 
holds for normal fluids, does not hold for polar and as- 
sociated fluids, has proposed recently to use the four pa- 
rameters T c , p c , Z c , and u>. This latter description can be 
then accounted for by using a four-parameter equation of 
state, but noting that the results obtained from the four- 
parameter equations which are constrained to reproduce 
the critical point, are only slightly better than those ob- 
tained from the three-parameter equations of state. 

A notable exception, recently proposed by Kiselev 
and Ely (author?) [24], is the empirical implementa- 
tion of the one-parameter (represented by the Ginzburg 
number Gi (author?) [25]) crossover description in a 
generalized corresponding-states model which use the 
(four-parameter) Patel-Teja equation of state = 

Poc 

f PT (^r, lo, ZqJ^J to calculate the classical behavior 
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of the Helmholtz free energy far away from the critical 

PT - 

point (with the condition Vo ^ c = Z^ C T < |). In that 
approach, the experimental value of the critical molar 
volume v c replaces ^ Tc as volume unit, leading to a 

calculated critical pressure p^J = Zq^^-^ from an em- 
pirical correlation expressing the calculated critical com- 
pression factor Z£ T (u, Z c ) as a function of real acentric 
factor ui and real critical compression factor Z c . A re- 
definition of the generalized corresponding-states princi- 
ple in the form -^r = / CR [Sr, ^-;u>, Z c ,Gi) only intro- 

Poc \ J o v " ) 

duces Gi as an additional corresponding-states param- 
eter. Hereabove / CR is a unique function, which ac- 
counts for a phenomenological crossover model that in- 
corporates singular behavior in the critical region, and 
transforms into an analytical equation of state far away 
from the critical point. Thus, it was assumed that the 
Ginzburg number Gi (uj, Z c , M mo i) can also be expressed 
as a function of real acentric factor u>, critical compres- 
sion factor Z c , and molar mass M mo \ of the fluid. As a fi- 
nal result, the unique four-parameter crosssover equation 

^It = / CR ffr-) Z^j, similar to the classical four- 
parameter equation of state, is able to predict with ac- 
ceptable accuracy the phase surface near to and far from 
the critical point, only using the real critical parameters 
T c , v c , Z C7 and the acentric factor lo. However, in such 

never be accounted for as entry parameters which char- 
acterize the one-component fluid. 

/From this brief status of the extended corresponding- 
states principle, it seems undeniable that a minimum 
set made of four parameters is necessary to character- 
ize each one-component fluid, thus identifying among a 
large set of uscfull fluid-dependent parameters the four 
critical parameters T c , p c , Z c , and <x c ,r (or Y c ). If the 
introduction of T c , p Cl v c , seems the natural way to de- 
fine the energy and length units and a related character- 
istic critical compression factor, the arbitrary choice of 
the fourth parameter confers empirical nature to any ex- 
tended corresponding-states approach (and to any func- 
tional form of the equation of state based on, especially 
when the calculated compression factor differs from the 
real one). 

It is then remarkable that the minimal set made of 
the four critical parameters T c , p c , Z Cl and Y c could be 
alternatively identified from a phenomenological analy- 
sis (author?) [26, 27, 28, 29, 30] of the singular be- 
haviour of these fluids approaching their liquid-gas crit- 
ical point. Moreoever, a fundamental distinction occurs 
in this critical phenomena scheme since the introduc- 
tion of the two dimcnsionless critical numbers Z c and 
Yc underlines their respective asymptotic scale-factor na- 
ture (see below and Refs. (author?) [26, 27]). As an 
essential consequence, the two-scale factor universality 
(author?) [31] estimated by the renormalization group 
method of field theory (author?) [32] is thus accounted 
for in building the master (i.e., unique) dimcnsionless 



crossover functions (author?) [33, 34, 35] of the renor- 
malized (i.e., rescaled) dimcnsionless field variables (au- 
thor?) [36, 37]. 

How to account for all these (practical and theoreti- 
cal) results to develop the appropriate forms of the an- 
cillary equations remains a difficult task, which was the 
object of few basic studies in regards the practical im- 
portance of the two-phase fluid properties. For example, 
the use of parachor correlations between a and Api,v 
are a convenient way of estimating surface tension from 
density measurements. Moreover, recent model calcula- 
tions (author?) [38] and phenomenological estimations 
(author?) [39] have shown the strength of the parachor 
correlations in the case of fluid mixtures. 

For pure fluids, parachor correlations have the follow- 
ing form 

where M mo \ = N Atrip. Na is Avogadro number and mp 
is the molecular mass. The subscript p refers to a molec- 
ular property, i.e. a property of the constitutive particle 
(atoms or molecules). The amplitude P e , called para- 
chor, is a fluid-dependent property, while the parachor 
exponent ir a ^ e is expected to have a unique numerical 
value for all the fluids. For a review of the n a ^ e and 
P e values, see for example Ref. (author?) [5]. A sub- 
script e indicates an effective value which corresponds to 
a given finite experimental range along the vapor-liquid 
equilibrium (VLE) line. Indeed, the introduction in the 
early 1920's of Eq. (1) was based on experimental ob- 
servations close to the triple point where pv ~ <C pl 
(and then Api,v ~ Pl), leading to the first proposed 
value of TT a .e = 4 for the effective parachor exponent 
(author?) [4]. It was then noted that the fluid para- 
chor P e is approximately a constant value in a large (in 
absolute scale) temperature range, leading to various at- 
tempts for its estimation by methods issued from group 
contribution methods, or extended corresponding-states 
principle. Accordingly, the parachor must be (at least) 
related to the four parameters involved in engineering 
equations of state. 

However, it is now well-established (author?) [40] 
that the validity range of such scaling form (1) is strictly 
restricted to the asymptotic approach of the liquid-gas 
critical point, where a and Ap^y simultaneously go to 
zero (author?) [41, 42] with the universal features of 
the uniaxial 3D Ising-like systems (author?) [31, 32]. 
Indeed, in that liquid-gas critical domain, it is expected 
that ApLy and a behave as Ap^y ~ (T c — T)' 3 and 
a ~ (T c —T) , respectively, where the critical exponents 
j3 and <j> take the following universal values at d = 3: 
~ 0.326 (author?) [43] and <j> = (d - 1) v * 1.26 (au- 
thor?) [41, 43]. To obtain the latter scaling law, use 
has been made of a ~ £ 1 (author?) [41, 42], where 
the power law £ ~ (T c - T)~ v (with v « 0.630 (au- 
thor?) [43]) accounts for the asymptotic singular bchav- 
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ior of the correlation length £. Obviously, the asymp- 
totic Ising-like form of Eq. (1) reads a oc (Ap^y) 3 , with 
To.e — > = |j ^ 3.87 when g — > and A/j^y — > 0, 
i.e., a universal power law behavior with universal criti- 
cal exponent 7r a of value which differs significantly from 
i"o,e = 4 for example. 

Therefore, despite the fact that the parachor corre- 
lations were initially developed to correlate the surface 
tension of liquids close to their triple point, their theo- 
retical justification is only well-understood close to the 
critical point. That infers a paradoxical situation when 
the main objective is to estimate a single value of P e in 
the largest two-phase domain. Moreover, when we con- 
sider an intermediate nonhomogeneous domain between 
critical point and triple point, it also seems not easy to 
invoke crossover arguments related to the classical mean- 
field theory of critical phenomena, as we will discuss in 
the final part of this paper. For example, effective values 
e ~ 0.36-0.30, (j) e ~ 1.20-1.30, and 7r a , e ~ 3.5-4, are 
observed in an extended temperature range which goes 
to the triple point (see Appendix A). Such values signif- 
icantly depart from mean-field ones /3mf = \, 4>mf = |, 
and 7r a: MF = 3 (author?) [44]. In addition to these dif- 
ferences, the mean-field exponents do not satisfy hyper- 
scaling (i.e. explicit d-dependence of some scaling laws), 
a difficulty precisely enhanced in the case of mean field 
exponents for the intcrfacial properties where an explicit 
(d — Independence also appears in "mixed" hyperscaling 
laws [for example, see below Eq. (51)]. 

Today, the estimation of parachors valid in a wide tem- 
perature range, from a limited number of fluid-dependent 
parameters, still seems an unsolved complex challenge. 

The main concern of the present paper is to clarify 
this situation by only using the four well-defined criti- 
cal parameters T c , p c , Z c and Y c in an asymptotic anal- 
ysis of the parachor correlations approaching the criti- 
cal point. A joint objective is to suppress asymptotical 
requirement for any other unknown adjustable param- 
eter in a well-defined extension of the VLE line close 
to the critical point, substituting then the two-scalc- 
factor universality of dimensonlcss fluids to the four- 
parameter corresponding-states principle to justify the 
observed master (i;e. unique) parachor function. 

The paper is organized as follows. In Section 2, we 
observe, by application of the scale dilatation method, 
the asymptotic master critical behaviors for interfacial 
properties close to the critical point. The scale dilata- 
tion method only uses a minimal set made of four criti- 
cal scale factors (neglecting here quantum effects in light 
fluids such as helium 3 (author?) [30]), noted Q" lm = 
j/?" 1 , a c , i? c , Y c } and defined in the next paragraph [see 
Eqs. (16) to (19)]. In Section 3, we unambiguously sepa- 
rate each scale factor contribution in the estimation of ei- 
ther the surface tension amplitude (only Y c -dependent), 
or the Ising-like parachor (only Z c -dependent). We con- 
clude in Section 4. Appendix A gives a complementary 
practical route to extend the analysis over the complete 
VLE line. 



2.1. 



MASTER SINGULAR BEHAVIORS OF 
INTERFACIAL PROPERTIES 

Asymptotic singular behavior of interfacial 
properties 



Close to the gas-liquid critical point, the asymptotic 
singular behavior of thermophysical properties are gener- 
ally characterized by Wcgncr-likc expansions (author?) 
[45] in terms of the following two relevant physical vari- 
ables (author?) [41, 46] 



At* = 



T-T c 



and 



P~ Pc 
Pc 



(2) 



(3) 



where subscript c refers to a critical property. At* and 
Ap are the temperature field and the order parameter 
density, respectively, of the liquid-gas transition. There- 
fore, in the non- homogeneous phase (At* < 0), along the 
critical isochore (p = p c ), the asymptotic singular behav- 
iors of the symmetrized order parameter density Ap^v 
and the intcrfacial tension a read as follows 



Ap LV = 2 Pc B | Ar* 



a = o- \At* 



1 



i=l 



|Ar 



* iiA 



(4) 



(5) 



The Ising-like universal values of the critical exponents 
are (3 ~ 0.326, <j> « 1.26, while A « 0.51 (author?) [43] 
is the Ising-like universal value of the lowest confluent 
exponent. The amplitudes B, Bi, cto, and <Xj, are fluid- 
dependent quantities which benefit from accurate the- 
oretical predictions of their universal combinations and 
universal ratios only valid in the Ising-like preasymtotic 
domain (author?) [47, 48], where the Wegner expan- 
sion are restricted to the first order term of the confluent 
corrections to scaling governed by the lowest confluent 
exponent A. 

Now, we introduce the squared capillary length (t-ca) 2 
[also called the Sugden factor (author?) [49, 50], noted 
S g ], related to a and Ap^y by 

2a 



V -Ca 



Sn 



g&PLv ' 



(6) 



where g is the gravitational acceleration. Equation (6) 
expresses the balance between interfacial forces and vo- 
lumic forces which define the shape and position of the 
liquid-gas interface in a gravity field of constant acceler- 
ation g. The asymptotic singular behavior of S g can be 
read as a Wegner-like expansion 



Sq = So \At* 



1 



_* iiA 



(7) 
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Figure 1: (Color online) a) Singular behavior (log-log scale) of 
the interfacial tension a (expressed in Nm" 1 ) and b) singular 
behavior of the symmetrized order parameter density ApLv 
(expressed in kgm~ 3 ), as a function of the temperature dis- 
tance T c — T (expressed in K) , for one-component fluids (see 
text, Appendix A and references in Table I). Inserted Table 
gives colour indexation of each selected fluid. The temper- 
ature axis is labelled by vertical arrows which indicate the 
practical values T = 0.99 T c and T = 0.7 T c , respectively (see 
text for details). 



leading to the scaling law 

<p = <f>-p 



(8) 



with if ~ 0.934, and to the canonical amplitude combi- 
nation 



So 



gpcB 



(9) 



Equations (4) to (9) are of basic interest as well for mea- 
surement techniques of interfacial properties as for theo- 
retical crossover descriptions when the temperature dis- 
tance to T c takes a finite value (see Appendix A). 

The singular behavior of the Sugden factor S g as a 
function of the temperature distance T c — T was illus- 
trated in Fig. 1 of Ref. (author?) [51] for about twenty 
pure compounds selected among inert gases, normal com- 
pounds and highly associating polar fluids. For the ref- 
erences of the S g and a data sources see the reviews of 
Rcfs. (author?) [51, 52, 53, 55]. Here, we have added 
the data sources (author?) [56, 57, 58, 59, 60, 61] of 



some hydrofluorocarbons (HFCs) and hydrocholorofluo- 
rocarbons (HCFCs) for related discussion in Appendix 
A. The raw data for the surface tension a (T c — T) and 
the symmetrized order parameter density Ap^y (T c — T) 
are reported in Figs, la and lb, respectively. For the 
ApLv data sources see for example the references given 
in Rcfs. (author?) [5], (author?) [55], and (author?) 
[56, 57, 58, 59, 60, 61]. In each case, the universal 
Ising-like slope of the asymptotic singular behavior ap- 
pears compatible with the experimental results. From 
these figures, it is also expected about one-decade varia- 
tion for each fluid-dependent amplitude 2p c B (T c ) and 
(To (T c ) ^ of the leading terms of Eqs. (4) and (5), re- 
spectively. 

As in the Sugden factor case (author?) [51], it also 
appears evident that the raw data used in present Fig. 1 
cover a large temperature range of the coexisting liquid- 
vapor phases since this range approaches the triple point 
temperature Trp . We have then adopted the same prac- 
tical distinction between asymptotic critical range and 
triple point region in the temperature axis, using vertical 
arrows for the temperature distances where T = 0.7 T c 
(i.e. the temperature distance where the fluid-dependent 
acentric factor u> is defined) and T = 0.99 T c . In a large 
temperature range defined by 0.3 < |At*| < 1 — ^p 2 -, the 
nonuniversal nature of each fluid is certainly dominant 
(see for example Appendix A), while, in the temperature 
range |At*| < 0.01, the singular behavior descriptions by 
Wegncr like expansions, and asymptotic two-scale-factor 
universality of their restricted two term form, hold. 

Selecting a and Ap^v measurements at identical val- 
ues of T c — T, we have constructed the corresponding 
a-ApLv data pairs. The singular behavior of a (Ap^y) 
is illustrated in Fig. 2a, while the corresponding behav- 
ior of a (ApLv)~™ a as a function of Ap^y is given in 
Fig. 2b, as usually made to enlighten the contribution 
of the confluent corrections to scaling and to have bet- 
ter estimation of the uncertainty attached to the value 
of the leading amplitude. Simultaneously, from xenon to 
n-octane, we also underline about a three-decade varia- 
tion for the fluid-dependent amplitude of this quantity, 
leading to about a half-decade variation of the effective 
parachor (expressed in [j mole -1 m 37ra ~ 2 ] unit when a, 
Ap LVl and M mol are expressed in [Jm -2 ] (or [Nm -1 ] ), 

[kgm -3 ], and [kgmole -1 ], respectively). Accordingly, 
by straightforward elimination of |At*| in Eqs. (4) to 
(7), we find the exact Ising-like asymptotic form 



{A PLV )'^Y 



(10) 



l + P l {A PLV y +0 
where (as previously mentioned in the introduction part) 

(11) 



^ = 1 = ^ + 1*3.87, 



The Ising-like asymptotic value Pq of the effective para- 
chor P e can be estimated using the leading terms of Eqs. 
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Figure 2: (Color on line) a) Asymptotic singular behavior of 
a as a function of ApLv obtained from Fig. 1 b) As in a) for 
the confluent quantity r^pfyW^ ( see ca Pti° n of Fig- 1 an d 
text for details and inserted Table for fluid colour indexation). 



the needed critical information to calculate any nonuni- 
versal leading amplitude of the selected fluid (here ne- 
glecting the quantum effects (author?) [30] to simplify 
the presentation) . The mass m,p of each molecule is also 
hypothesized known to infer the total amount N of fluid 
particles by measurements of the fluid total mass M = 
NAmp. p (p c ) is the (critical) pressure. Vp = -^j = — 

(vp.c = is the molecular volume (critical volume). 

Nvp is the extensive variable con- 



The total volume V 
jugated to p. 7 C = 



(0p\ 
\dT J 



is the common critical 



CP 



direction in the p; T diagram of the critical isochore and 
the saturation pressure curve p sat (T) at critical point 
(CP), thus defined by 



7 C = 



dp 
df 



dp s 



Rewriting Eq. (13) as a four-scale-factor set 
Qf n = {f3- 1 ,a c ,Z c ,Y c } 

where, 

Pc 1 = k B T c ~ [energy] 



[length 



(14) 

(15) 
(16) 
(17) 



(4) to (7) and reads 



, .9 M mol 



At, 



mol 



2{p c B) 1 * 



(12) 



depending on the pair of selected variables, either 
{a; ApLv} or {S g ; ApLv}- Equations (10) and (11) 
clearly demonstrate the critical scaling nature of Eq. 
(1), with an essential consequence: the parachor Pq is 
a nonuniversal leading amplitude which must satisfy the 
two-scale-factor universality of the Ising-like universality 
class. 

How to estimate the parachor appears thus as a ba- 
sic question in a sense that only two nonuniversal lead- 
ing amplitudes arc sufficient to characterize the complete 
singular behavior of a one-component fluid when T — > T c 
and Ap^ — > 0. 



2.2. The basic set of fluid-dependent parameters 

As proposed by Garrabos (author?) [26, 27], a phe- 
nomenological response to the above basic question relies 
on the hypothesis that the set 



x-={pc,%,c,t C)7 :}, 



(13) 



of four critical coordinates which localize the gas-liquid 
critical point on the p, i%, T phase surface, contains all 



Z r = 



p c m p 
ksTcPc 



Yc = l c - 

Pc 



(18) 



(19) 



we introduce the energy unit [/3c -1 ], the length unit [a c ], 
the (isothermal) scale factor [Z c ] of the order parameter 
density [see Eq. (3)] along the critical isothermal line, 
and the (isochoric) scale factor [Y c ] of the temperature 
field [see Eq. (2)] along the critical isochoric line. 

Table I provides values of the critical parameters in- 
volved in Eqs. (13) and (15) for the 23 pure fluids se- 
lected in this paper. 

Indeed, our dimensional scale units of energy (/3 C ) _1 
[Eq. (16)], and length a c [Eq. (17)], provide equiva- 
lent description as T c and p c of the basic (2-parameter) 
corresponding-states principle. The customary dimcn- 
sionless forms of a and Ap^y are a* = (a c ) d_1 /3 c a and 
ApLv = A 2p V ■ A s for the Sugden factor case, Fig. 

3 gives distinct curves of a* (ApLv) ™ a as a function 
of Ap~Lv, which confirms the failure of any description 
based on the two-parameter corresponding-states prin- 
ciple. Moreover, the direction difference 7r aj MF - fa ^ 
—0.87 with a classical power law of mean field exponent 
7i"a,MF = 3 also disagrees with experimental trends at 
large temperature distance, as illustrated in Fig. 3. 
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Fluid 


mp 


T c 


Pc 


pc 


7c 


(/3c) 1 


Ctc 


z c 


Y c 




(io- 2e k g ) 


(K) 


(MPa) 


(kgm- 3 ) 


(MPaK" 1 ) 


(nr 21 j) 


H 






Ar 


6.6335 


150.725 


4.865 


535 


0.19025 


2.08099 


0.7535 


0.289871 


4.89423 


Kr 


13.9153 


209.286 


5.500 


910 


0.1562 


2.8895 


0.8069 


0.291065 


4.94372 


Xe 


21.8050 


289.740 


5.840 


1113 


0.1185 


4.0003 


0.8815 


0.286010 


4.87914 


N 2 


4.6517 


126.200 


3.400 


314 


0.1645 


1.74258 


0.8002 


0.289078 


5.10585 


o 2 


5.3136 


154.58 


5.043 


436 


0.1953 


2.13421 


0.7508 


0.287972 


4.98641 


C0 2 


7.3080 


304.107 


7.732 


467.8 


0.173 


4.19907 


0.8289 


0.274352 


6.13653 


SF 6 


24.2555 


318.687 


3.76 


741.5 


0.084 


4.40062 


1.054 


0.281243 


6.11960 


H 2 


2.9969 


647.067 


22.046 


322.8 


0.2676319 


8.93373 


0.7400 


0.229117 


6.85520 


C2H4 


4.658 


282.345 


5.042 


214.5 


0.11337 


3.89820 


0.91781 


0.28131 


5.34856 


CH 4 


2.6640 


190.564 


4.59920 


162.7 


0.14746 


2.63102 


0.8301 


0.285752 


4.981927 


C2H6 


4.99324 


305.322 


4.872 


206.58 


0.10304 


4.21554 


0.95290 


0.27935 


5.45505 


C3H8 


7.32248 


369.825 


4.2462 


220 


0.077006 


5.106 


1.063 


0.27679 


5.70688 


11-C4H10 


9.6518 


425.38 


3.809 


229 


0.0643 


5.87301 


1.155 


0.273352 


6.17774 


i-C4Hio 


9.6518 


407.85 


3.65 


225 


0.0643 


5.63102 


1.155 


0.278056 


6.18173 


C5H12 


11.9808 


469.70 


3.3665 


232 


0.0511 


6.48491 


1.244 


0.270875 


6.12956 


CeHi4 


14.3100 


507.85 


3.0181 


234 


0.043658 


7.00666 


1.319 


0.266670 


6.30719 


C7H16 


16.6393 


540.13 


2.727 


234 


0.038068 


7.45731 


1.398 


0.262180 


6.64356 


CgHig 


18.9685 


568.88 


2.486 


232 


0.033768 


7.85424 


1.467 


0.258978 


6.82776 


HFC-32 


8.6386 


351.26 


5.782 


423 


0.124088 


4.849676 


0.9431 


0.243384 


6.53842 


HCFC-123 


25.3948 


456.82 


3.666 


554 


0.005711 


6.30764 


1.1982 


0.266439 


6.11647 


HCFC-124 


22.6622 


395.35 


3.615 


566 


0.063205 


5.45840 


1.1473 


0.265179 


5.91234 


HFC- 125 


19.9301 


339.17 


3.618 


568 


0.078665 


4.68275 


1.0898 


0.271100 


6.37446 


HFC-134a 


16.9426 


374.30 


4.065 


512.7 


0.083109 


5.16777 


1.0833 


0.259940 


6.65260 


HCFC-141b 


19.42 


477.31 


4.250 


460 


0.060522 


6.58998 


1.1575 


0.272268 


5.79119 


HCFC-142b 


16.6876 


410.26 


4.041 


447 


0.071972 


5.66426 


1.1191 


0.266338 


6.30691 


HFC-152a 


10.9680 


386.41 


4.512 


369 


0.086345 


5.33497 


1.0574 


0.251384 


6.39463 



Table I: Critical parameters for the selected one-components fluids. 
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Figure 3: (Color online) Singular behavior (log-log scale) of 
the dimensionless quantity ^-^ j,, as a function of the sym- 
metrized order parameter density A/5_lv; inserted slope: di- 
rection difference 7T ,mf — 7r a = —0.87 from 7t q ,mf = 3 (see 
text); arrow sets in lower horizontal axis: I At* I = 10~ 2 and 
At* I = 0.3 (see text and caption of 1); Inserted Table gives 
the fluid colour indexation. 



On the other hand, the set Q™ n [Eqs. (15) to (19)] 
conforms to the general description provided from ap- 
propriate 4-parameter corresponding-states modelling, as 
mentioned in our introductive part. For example, we re- 
trieve that the scale factor Y c [Eq. (19)] is related to the 
Riedel factor a c .p — 7 C ^ by Y c + 1 = a c .n, as mentioned 
in our introductive part. However, as we will extensively 
show in this paper, the scale dilatation approach brings a 
theoretical justification to these critical parameters, ini- 
tially introduced only to build extended corresponding- 
states principle. 

First, the microscopic meaning of (/3 C ) and a c , re- 
lated to the minimum value [e~(/? c ) _1 ], of the interac- 
tion energy between particle pairs at equilibrium position 
[r^iaj , takes primary importance. Thus, a c appears 
as the mean value of the finite range of the attractive 
interaction forces between particles. So that, 



v c ,i 



(20) 



is the microscopic volume of the critical interaction cell 
at the exact (T = T c and p — p c ) critical point. Us- 
ing thermodynamic properties per particle, it is straight 



9 



to show that the critical interaction cell is filled by the 
critical number of particles 



1 



since, rewriting Eq. (18), we have 



Pc 



V c ,i = N cJ x v Pt> 



X — - 

Z c Pc 



(21) 



(22) 



Second, since in the critical phenomena description only 
one length scale unit is needed to express correctly the 
non-trivial length dimensions of thermodynamic and cor- 
relation variables (author?) [31], the above microscopic 
analysis is of primary importance. By chosing a c , the 
two dimensionless scale factors Z c and Y c are then char- 
acteristic properties of the critical interaction cell of each 
one-component fluid. Specially, -A- takes similar micro- 
scopic nature of the coordination number in the lattice 
description of the three-dimensional Ising systems, while 
a c takes similar microscopic nature of their lattice spac- 
ing dlsing. 

On the basis of this microscopic understanding, we are 
now in position to estimate the Ising-likc parachor from 
<2™ m , using scale dilatation of the physical fields (au- 
thor?) [26]. 



2.3. The scale dilatation of the physical variables 

The asymptotic master critical behavior for interfacial 
properties when T — > T c and Ap^y — > 0, can be observed 
by using the following dimensionless physical quantities 

At* = k B (3 c (T c - T) (23) 
Ap* = p c (fi p - n P , c ) (24) 

Am* = (a c ) d (n~n c ) = (Z^ 1 Ap (25) 
and the following master (rescaled) quantities, 

T* =Y c \At*\ (26) 



g* = rripP c a c g. pp (pp, c ) is the molecular chemical po- 
tential (critical molecular chemical potential), n (n c ) is 
the number density (critical number density). The nor- 
malized variable n = where pp is conjugated to the 
total amount of matter N, is then related to the order 
parameter (number) density expressing thermodynamic 
properties per molecule. As mentionned in the introduc- 
tion, the (mass) density p = where p p is conjugated 
to the total mass of matter M, is related to the order 
parameter density, but expressing thermodynamic prop- 
erties per volume unit. Here we have noted p p = -^f- 
the chemical potential per mass unit. The dimension- 



less form of p p is p p 



while the one of 



Pp is 



have used a c (mpf3 c )~- 



as a time unit. 



Pp = ppftc, with -g-Pp = p P - g is the gravitational accel- 
eration needed to perform conventional measurements by 
capillary rise or drop techniques (author?) [2, 49]. g* 
is the dimensionless gravitational acceleration, where we 

mass(lcngth) 2 
{ energy 

In Eqs. (26) and (27), Y c and Z c are two independent 
scale factors that dilate the temperature field along the 
critical isochore and the ordering field along the critical 
isotherm, respectively. These Eqs. (26) and (27) are for- 
mally analogous to analytical relations (author?) [62] 
linking two relevant fields of the so-called $j_ 3 (n = 1) 
model with two physical variables of a real system belong- 
ing to the universality class. That implicitcly imposes 
that only a single microscopic length is characteristic of 
the system (author?) [31], which then can be related 
to the inverse coupling constant of the model taking ap- 
propriate length dimension, precisely for the d — 3 case. 
When a single length is the common unit to the thermo- 
dynamic and correlation functions, the singular part of 
free energy of any system belonging to the universality 
class remains proportional to a universal quantity which 
generally refers to the (physical) value of the critical tem- 
perature. In the case of one-component fluids, the sin- 
gle characteristic length originates from thermodynamic 
considerations [see above a c of Eq. (17)]. Accordingly, 
the universal singular free energy density is expressed in 
units of (/3 C ) _1 = ksT c [sec Eq. (16)]. The master fields 
T* and TL* have similar Ising-like nature to the renor- 
malized fields t and h in field theory applied to critical 
phenomena. 



U* = {Z c y- Ap P = {Z c )-- A~p p 



M* = {Z c y Am* = (Z c ) 5 Ap 



S*=g* {Z c )-*{t Ca ) 



d-l 



(27) 



(28) 



(29) 



(30) 



where, in Eqs. (26) to (30), we have only used the two 
scale factors Y c and Z c to rescale the dimensionless quan- 
tities, with (j* = (ct c ) d 1 ftcO", P*ca = ( a c) 1 ^Ca and 



2.4. Master crossover behavior for interfacial 
properties of the one-component fluid subclass 

The observation of the critical crossover behavior of 
a master property V* in a V* — T* diagram, gener- 
ates a single curve which can be described by a mas- 
ter Wegner-like expansion V* (T*). Asymptotically, i.e. 
for T* — » 0, the universal features of fluid singular be- 
haviors are at least valid in the Ising-like preasymptotic 
domain (author?) [34, 47] where only two asymptotic 
amplitudes and one first confluent amplitude character- 
ize each one-component fluid. In our present formu- 
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Fluid 


(f> e 


O"0,e 


Ref. 


00 


2><T, CXp 




Ref 






(lO^Nm" 1 ) 




(lO^Nm- 1 ) 




% 




Ar 








31.55 


1.162 


0.5 


(author?) [60] 










30.48 


1.123 


-2.9 


(author?) [52] 




1.281 


38.07 


(author?) [63] 


31.42 


1.158 


0.1 


This work 




1.277 


37.78 


(author?) [65] 


31.76 


1.170 


1.2 


This work 


Xe 








45.69 


1.179 


2.0 


(author?) [53] 










45.61 


1.177 


1.8 


(author?) [52] 




1.302 (±0.006) 


62.9 (±1.8) 


(author?) [67] 


47.1 


1.215 


5.1 


This work 




1.290 


53.9 


(author?) [66] (from (author?) [64]) 


42.68 


1.101 


-4.8 


This work 




1.287 (±0.017) 


54.6 (±0.1) 


(author?) [66] 


43.83 


1.131 


-2.2 


This work 


N 2 








25.75 


1.212 


4.8 


(author?) [52] 


2 








34.04 


1.186 


2.6 


(author?) [52] 


C0 2 








70.05 


1.209 


4.6 


(author?) [53] 










65.93 


1.138 


-1.5 


(author?) [52] 




1.26 


76 


(author?) [52] (from (author?) [68]) 


69.1 


1.193 


3.2 


This work 




1.281 


84.72 


(author?) [70] 


69.92 


1.207 


4.4 




SF 6 








47.85 


1.230 


6.4 


(author?) [53] 










44.09 


1.134 


-1.9 


(author?) [52] 




1.285 (±0;016) 


55.13 (±2.6) 


(author?) [69] 


46 


1.183 


2.3 


(author?) [52] 




1.285 (±0;016) 


55.13 (±2.6) 


(author?) [69] 


44.64 


1.148 


-0.7 


This work 




1.286 


54.28 


(author?) [71] 


43.78 


1.126 


-2.6 


This work 


CBrF 3 


1.279 


54.05 


(author?) [71] 


45.02 


1.106 


-4.3 




CCIF3 


1.30 


58.84 


(author?) [70] 


44.5 


1.135 


-1.8 


This work 




1.283 


52.53 


(author?) [71] 


42.96 


1.096 


-5.2 


This work 


CHCIF2 


1.283 


69.03 


(author?) [71] 


56.44 


1.139 


-1.5 


This work 


CCI2F2 


1.283 


59.63 


(author?) [71] 


48.76 


1.172 


1.4 


This work 


CCI3F 


1.263 


63.24 


(author?) [71] 


56.78 


1.196 


3.5 


This work 


H 2 








218 


1.135 


-1.8 


(author?) [53] 










220.7 


1.149 


-0.6 


(author?) [52] 


CH4 








30.22 


1.012 


-12.4 


(author?) [52] 


C2H4 








45.08 


1.166 


0.8 


(author?) [53] 


C2H6 








45.15 


1.145 


-1.0 


(author?) [53] 


i-C4Hio 








45.9 


1.159 


0.2 


(author?) [60] 


mean 










1.167 


1.0 




standard 










0.062 


5.4 





Table II: Effective values of the critical exponent (column 2) and asymptotic amplitudes of interfacial tension (column 3) from 
references given in column 4. Related estimations of the leading asymptotic amplitude 00 (column 5) (see text for details). 
Calculated values of the corresponding master amplitudes -Z CT ,ex P = <jq (ct c ) 2 (3 C (Yc) - '* (column 6) of the master interfacial 
tension (see Table I for the values of a c , (3 C , and Y c ). The % differences 8Z a . c ^ p — 100 ^ z y*p _ )J f rom the values Z a = 1.156 
estimated from universal amplitude combinations are given in column 7. For references see column 8. 



lation of the liquid-vapor interfacial properties, we de- 



LV 



(Z c y Am* LV = (Z c )*Ap LV , Am 



LV 



2/', 



and 



fine M 

(a c ) d (n L -n v ) = (Z C )~ L Ap LV , Ap LV 
we introduce the following master equations of interest, 
here restricted to the first order of the critical confluent 
correction to scaling to be in conformity with the above 
universal facturcs, 



M* LV = Z M {T*Y 1 + Z$(T*) 



(31) 



E* = Z a {T*) 9 1±Z«(T*) A 



s;, =z s {T*y fi + zW (T*) 

with their interrelation [see Eqs. (8) and (6)] 
9* - S * 

JVI LV 



(32) 
(33) 

(34) 
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The master asymptotic behaviors of M* LV and S*» as 
a function of T* were observed and analyzed in Refs. 
(author?) [29] and (author?) [51] for several pure flu- 
ids. The corresponding master amplitudes take the val- 
ues Z M ~ 0.468 (±0.001) and Z s ~ 2.47 (±0.17) (for the 
quoted error-bars see Refs. (author?) [29, 51]). /From 
Eq. (34), the leading master amplitude for the surface 
tension case is then Z a = ZmZs ^ 1.156 (±0.087). 
Here, using a method similar to the one applied to the 
Sugden factor case, we show in Table II that this mas- 
ter value is compatible with the results obtained from 
experiments performed sufficiently close to the critical 
point. Therefore, only surface tension measurements al- 
ready analyzed by Moldover (author?) [53] and Gie- 
len et al (author?) [52] are considered. The respective 
effective values 4> e and <ro,e of the exponent-amplitude 
pair corresponding to data fits by an effective power law 
o = co,e IAt*!^" are reported in columns 2 and 3, when 
needed for the present analysis (see the corresponding 
Refs. (author?) [63, 64, 65, 66, 67, 68, 69, 70, 71] in 
column 4). The estimated values of the leading ampli- 
tude (To for the Ising value <j> ~ 1.260 of the critical expo- 
nent are given in column 5. The references reported in 
column 8 precise the origin of these estimations, which 
are mainly dependent on the accuracy of the interfacial 
property measurements in the vicinity of |Ar*| ~ 0.01. 
For example, in the present work we have estimated 

0o by the following relation do = 



l.i 



us- 



ing data sources of columns 2 and 3. Our estimation is 
then compatible with surface tension measurements at 
At* | = 0.01, neglecting confluent corrections in Sug- 
den factor measurements (author?) [51], and averag- 
ing (for all the selected fluids) the confluent correction 
contributions in density measurements to 10% at this fi- 
nite distance to the critical temperature (author?) [29]. 
The values of Z a ^ v = <jq (a c ) 2 I3 C {Y C )-* (column 6) 
calculated from these Co estimations arc in close agree- 
ment with our master value Z a = 1.156. The residuals 

SZ aiCxp = 100 (^f 31 - l) (column 7) expressed in %, 

are of the same order of magnitude than the experimen- 
tal uncertainties [see for example Refs. (author?) [53] 
and (author?) [52]]. We note that the ±1% residuals 
between the experimental mean value {Z ajfSX p) = 1.167 
and the estimated master one Z a = 1.156, have a stan- 
dard deviation (±5.4%) comparable to the experimental 
uncertainties (~ 7.5%). However, this good agreement 
on the central value is noticeable in regards to signifi- 
cant contributions of confluent corrections, as reflected 
by the effective exponent values 4> e « 1.28 — 1.31 larger 
than (f> ~ 1.260 at finite distance from T c (see also a com- 
plementary discussion related to the analysis of data at 
large distance from T c given in Appendix A). 



Similarly, each confluent amplitude Z 

r(l) 



(1) 
M ' 



and 



Despite their great interest for the validation of the- 
oretical predictions, the exact forms of Eqs. (31) to 
(33) are not essential to understand the independent 
scaling role of the two scale factors Y c and Z c . More- 
over, our present interest is mainly focused on the scal- 
ing nature of Z c . More precisely, by a variable exchange 
i 

from A/5 to (Z c ) 2 Ap, the expected master behavior must 

be observed when the master order parameter density 

i 

M* = (Z c ) 2 Ap is used as a x- axis, generating then 
each "collapsed" curve of master equation V* (M*) in 
the Ising-like preasymptotic domain, i.e. for M* — > 0. 
Usually, such a master behavior of the singular fluid 
(bulk) properties occurs along the critical isothermal line 
T = T c , i.e., for T* = At* = 0, which provides disap- 
pearing of the scale factor Y c and which only preserves the 
contribution of the scale factor Z c in the determination 
of the fluid-dependent amplitudes (author?) [27, 28]. 
However, in the nonhomogeneous domain, the order pa- 
rameter density spontaneously takes a finite value to dis- 
tinguish the two coexisting phases in equilibrium. Thus, 
along the critical isochore, it is also possible to observe 
the master behavior of any singular (interfacial) prop- 
erty expressed as a function of the "symmetrized" order 
parameter density, starting with the surface tension as a 
typical example. 

/From a and Ap^y at identical T c — T, we can 
construct data points of master coordinates £* = 

(Z c )^ Ae^v i n the Y,*-M*, 



(a c ) d 1 Pc® and A4 



l LV ~ y^cj 2 p c 111 ^ J l LV 

diagram. As illustrated in Fig. 4a, all these data 
collapse to define a master behavior of £* (M* LV ) for 
M 



LV 



0. This collapse is well-enlightened in Fig. 4b 
which illustrates the corresponding master behavior of 
S* (M* LV )~ 7ra , without any reference to a fitting mas- 
ter equation. To have a better appreciation of the real 
temperature range of the VLE domain, the values of 
| At* | = 0.01 and |At*| = 0.3 arc also given by (fluid- 
dependent) arrows in the M* LV axis. 

Obviously, using Wegner-like expansions to eliminate 
T* between Eqs. (31) and (32), the leading power law of 
S* as a function of M* LV reads as follows 

S* = Z a (M* LV )% [l + o[{M* Lv y]) (35) 



where Z„ 



~ 21.88. This asymptotical am- 



plitudc is illustrated by an arrow in the vertical axis of 
Fig. 4b and the related horizontal dashed (blue) line 
indicates clearly that the observed master behavior at fi- 
nite distance to the critical point seems in "asymptotical" 
agreement. 

In addition, the formal analogy between the basic hy- 
potheses of the rcnormalization theory and the scale di- 
latation method makes also easy to probe that the effec- 
tive fluid crossover behavior estimated from the massive 
rcnormalization scheme is consistent with the asymptoti- 



Z y s , takes a master constant value for all one-component cal power law of Eq ^ ( 35 )] when m * l q. As a matter 



fluids. Among all Zf J , only one is independent and char- 
acteristic of the pure fluid subclass. 



LV 

of fact, along the critical isochore, the theoretical esti- 
mations of the fluid master behaviors in a whole thermal 
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Figure 4: (Colour online) a) Master singular behavior (log- log scale) of the renormalized surface tension E* = a* [Eq. (29)], 
as a function of the renormalized (symmetrized) order parameter density M* LV in the nonhomogeneous domain [see Eq. (28) 
and text]; b) As in a) for the "confluent" quantity -, — f (see text). In a) and b): Dashed blue curve: Eq. (36); full red 

[ M LVJ 

curve: (37); full arrow with label PAD:^^ < -^pad extension of the preasymptotic domain; dotted arrow with label EAD: 
M*lv < M^ 1 ^ extension of the extended asymptotic domain; The graduation of the upper horizontal axis gives i* (T* < 0) 
calculated from theoretical crossover (see text and Ref. (author?) [51]); arrow sets in lower horizontal axis: | At* | = 1CT 2 and 
| At* | = 0.3, respectively; Inserted Table gives the fluid color indexation. 



field range < |T*| < oo can be made using the recent 
modifications (author?) [34, 35] of the crossover func- 
tions calculated by Bagnuls and Bcrvillicr (author?) [33] 
for the classical-to-critical crossover of the Ising-like uni- 



versality class. However, as already noted in the Sug- 
den factor case (author?) [51], the theoretical function 
giving the classical-to-critical crossover of the interfacial 
tension is not available. Thus, to estimate E* h , we must 
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use two alternative routes, introducing the theoretical es- 
timations of either M* h (T* < 0) or C* h (T* < 0). 

A first straightforward route consists in using as 
an entry quantity the theoretical value A4* h of the 
master order parameter density, admitting then that 

Ml h (T < 0) = M* LV (|T|), an identity already validated 

■f l f \ 

at least in the preasymptotic domain |T| < ^pad ( see 
Ref. (author?) [29]). As a result, the theoretical para- 
chor function reads as follows 



J th,M 



Za (M* th ) 



(36) 



justifying precisely the dashed blue curves in Figs. 4a 
and b, for the complete M* LV -va,ngc. Unfortunately, at 
large values of T* , the theoretical function M* th (T* < 0) 
is certainly not able to reproduce the experimental be- 
havior of A^y(|T|) and, in the absence of this dedi- 
cated analysis, we cannot indicate easily the true exten- 
sion of the M* LV -T&nge where the identity M* h (T < 0) = 
M* LV (\T\) holds. 

To pass round this difficulty, a second route combines 
the theoretical functions £* h (T* < 0) and M* th (T* < 0) 
in a whole thermal field range < \T* \ < oo, to infer nu- 
merically the function £* h (Ai* h ) by exchanging the |T*|- 
dependcncc of £* h by the A^-depcndcncc of |T*| (i.e., 
reversing the function {T* < 0)). Anticipating the 
introduction of the universal number recalled below, 
the following asymptotic scaling form of the renormalized 
interfacial tension 



J thJ 



KeltUMtk)] 1 - 



(37) 



provides the correct asymptotic behavior of £* h (Ml h ) 
in the Ising limit A4* h — > 0, as illustrated by the full 
red curves in Figs. 4a and b. However, in spite of the 
increasing difference observed in Fig. 4 between these 
two theoretical estimations, or between them and the ob- 
served master behavior, when M* LV increase, we can now 
give a well-defined estimation of the effective extension 
of the VLE domain where the master crossover behavior 
of E*have physical meaning. 

Indeed, in a similar way as in our previous analy- 
sis of S g , the universal prefactor of Eq. (37) accounts 
for the "Ising-like" behaviors of the correlation length 
£ (At* < 0) and the surface tension a (|Ar*|), where we 
have used the universal ratio |[at*<o) = 1.96 and the 
universal amplitude combinations given by the prod- 
ucts of the interfacial tension by the squared correlation 
length (author?) [52, 53], i.e., 



R%=]im{p c o(\Ar*\)[t;(AT*)] d - 1 } ATt 



(38) 



where « 0.376 = (1.96) R~ v so that R~ £ « 0.979 
(for the estimated error-bars see also Ref. (author?) 
[51]). The superscripts ± refer to the singular behavior 
of £ above (+) or below (— ) T c . Equation (38) means 
that the interfacial energy of a surface area tends 
to a universal value [expressed in units of (/3 C ) ] for any 



system belonging to the Ising-like universality class (we 
recall that the thickness of the interface is then of or- 
der £ (At* < 0)). In such a description, the Wegner-like 
expansion of the correlation length, along the critical iso- 
chore, above and below T r , reads as follows 



1+X> I 



(39) 



where the fluid-dependent amplitudes £q and £ are such 

that fir ~ 1.96 (see above), while the contribution of 

the confluent corrections is hypothesized the same above 
and below T c . Accordingly, the master correlation length 
£* h = — (where we neglect here the quantum effects at 
the microscopic length scale of the order of a c (author?) 
[30]), can be estimated using the asymptotical modifica- 
tions (author?) [35, 72] of the crossover function for 
the correlation length in the homogeneous domain (au- 
thor?) [33]. This master asymptotic behavior of I* as 
a function of T* > was analyzed in Ref. (author?) 
[72] for seven different pure fluids, demonstrating that 
two specific £* h -values in the range £^ h > 1 provide con- 
venient marks to define: 

i) the extension of the Ising-like preasymptotic domain 
(i.e., £* h > 80), where the fluid characterization is ex- 
actly conform to the universal features calculated from 
the massive rcnormalization scheme; 

ii) the extension of the effective fluid master bahavior 
at finite distance to the critical temperature where £% h > 



_ r th (T*>o) 



we can then 



Introducing £* th (T* < 0) = 
complete the restricted two-term master forms of Eqs. 
(31) to (33) valid in the Ising-like preasymptotic domain, 
by the following two-term equation 



t (T* < 0) = Z~ \T*\ 



1 



Z\ \T* 



(40) 



where iJ, 



1.96 



j ~ 0.291 (with Z+ « 0.570) and 

Z\ = Z\' + = 0.377 (author?) [72]. The universal fea- 
tures of the interfacial properties within the Ising-like 
preasymptotic domain are in conformity with the three- 
master amplitude characterization of the one-component 
fluid subclass defined in Ref. (author?) [35]. The 
related singular behaviors of the dimensionlcss interfa- 
cial properties of each pure fluid can be estimated only 
knowing Y c and Z c . Using now the numerical function 
£j* h (A4* h ), we have also normed the upper x-axis in Fig. 
4 to illustrate the singular divergence of (T* < 0) in 
complete equivalence to the upper x-axis of Fig. 3 in Ref. 
(author?) [51]. 

iFrom |T*| = Z^ad - 5 x 10 ~ 4 where e th - 40 [ see 
Fig. 3 in Ref. (author?) [51]], the A^y-extension of 
the preasymptotic domain is 



M* LV < M§& 



4 x 10" 



(41) 



as shown by the arrow labeled PAD in Fig. 4. Within this 
Ising-like preasymptotic domain, the agreement between 
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the two theoretical estimations of S* arc (qualitatively) 
conform to the three-amplitude characterization of the 
universal features. The quantitative conformity cannot 
be exactly accounted for within a well-estimated theoret- 
ical erro-bar, due to the absence of theoretical prediction 
for the crossover of the surface tension, large error-bar in 
the estimation of the amplitude of the first-order conflu- 
ent correction term of the order-parameter density, and 
hypothesized contribution of the (homogeneous) conflu- 
ent corrections in the correlation length case. Since a c 
measures the shorted-range of the microscopic molecular 
interaction, £*=-£- gives the relative order of magni- 
tude of the true correlation length £, and we can retrieve 
for the properties of a vapor-liquid interface of thickness 
~ £, the similar physical meaning of the two-scale factor 
universality in the close vicinity of the critical point, i.e., 
when the conditions £ 3> a c , or equivalcntly £* 3> 1, 
are satisfied. Practically, the asymptotic singular be- 
haviors of a (respectively £*) and &.p* LV (respectively 
M* LV ), including then the first order confluent correc- 
tion to scaling as given by Eqs. (31) to (33), are observed 
when the correlation length in the non-homogeneous do- 
main estimated from Eq. (40) is such that £ > 40 a c , or 

{M* LV < 7\4ad - °- 04 ; l T *l £ 4ad - 0.0005}, equiv- 
alcntly. This finite extension of the non-homogeneous 
Ising-like preasymptotic domain, corresponds to a cor- 
relation volume £ d of the fluctuating interface (of typi- 
cal thickness ~ 40 nm, sec Tabic I) which contains more 
than 6x 10 4 "microscopic" (i.e., v c j) volumes, and there- 
fore, at least 2 x 10 5 cooperative particles for which the 
microscopic details of their molecular interaction at the 
a c -scale (typically ~ 1 nra, see I) are then unimportant 
(here admitting that mean number of fluid particles fill- 
ing v c j is typically -J- ~ 3.5). 

Simarly, using \T*\ = £W> ~ 0.03 [see Eq. (55) in 
Ref. (author?) [51]] where £j* h ~ 3, the M* LV -cxtensicm 
of the extended asymptotic domain is 

M* LV < M {lf} -0.16 (42) 

as shown by the arrow labeled EAD in Fig. 4). As ex- 
pected, the master behavior of E* is readily observed in 
this extended critical domain and in Appendix A, we give 
a convenient master modification of Eq. (36) to account 
for it precisely. Such an extended domain of the mas- 
ter behavior of the fluid subclass can be still understood 
since £ > 3a c , so that £ d > 30v c j, and then more than 
100 particles in cooperative interaction. However, it is 
noticeable that the practical relative values T = 0.99 T c 
of the temperature distance to T c frequently referred to 
define the critical region for each pure fluid are not inside 
the effective extension of the observed master singular 
behavior (see also Appendix A). 

Finally, using \T*\ = £cic — 0.2 where the size of 
the correlation length is equal to the size of the critical 
interaction cell filled by three or four particles, i.e., £* h ~ 
1 [see Fig. 3 in Ref. (author?) [51]], the related value 



of the order parameter density is 

M C ic - 0.3 (43) 

as illustrated by vertical dotted line in Fig. 4. We note 
that the pratical limit T = 0.99 T c with £ ^ 2a c , is 
in between MP~^ and Mcic, so that £ d ~ Sv c j only 
involving ~ 28 particles in interaction. Such a micro- 
scopic situation make questionable the critical nature of 
the fluid properties measured at this finite distance to 
T c , and more generally, shows that the range M* LV > 0.3 
can be considered as non-Ising-like in nature. Especially, 
Fig. 4 indicates unambiguously that the temperature 
T = 0.7 T c where the acentric factor is defined, leads 
to a correlation length smaller than the mean equilib- 
rium distance r e between two- interacting particles, since 
the value I* ~ | (see the hatched limit labeled m in 
upper x-axis of Fig. 4) corresponds approximatively to 
£ ^ \ol c ^ r e . Appendix A provides complementary 
analysis of this "nonuniversal" fluid crossover over the 
complete temperature range. 

The remaining correlative problem applying the scale 
dilatation method to liquid-vapor interfacial measure- 
ments, is to estimate the respective contribution of each 
scale factor Y c and Z c in the fluid-dependent amplitudes 
of the surface tension, expressed either as a function of 
T* , or as a function of M* LV . This problem is treated in 
the next section. 



3. INDEPENDENT SCALING ROLES OF THE 
TWO SCALE FACTORS 

3.1. The thermal field dependence characterized by 
the Y c scale factor 

By inverting Eqs. (26) to (30), we can easily recover 
the asymptotical form for the interfacial properties of 
Eqs. (4) to (7) from the asymptotic form of the master 
Eqs. (31) to (33). For example, the leading physical 
amplitudes B, <jq, and So, can be estimated from the 
following relations 

B = {Z B )-*{Y c f Z M (44) 
a o = (/3 c )- 1 (a c ) 1 - d (r c ) Z CT (45) 

So = {a.f- 1 (g*)- 1 (Zcf 2 {Y c y Z s (46) 

Similarly, from comparison of Eqs. (39) and (40), the 
amplitudes ^ of the (bulk) correlation length can be 
estimated from the following relation 

it = a c {Y c )- v Zf (47) 

As expected, the leading amplitudes arc combinations of 
the fluid scale factors. However, we underline the uni- 
versal hyperscaling feature of Eqs. (45) and (47), where 
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<7o and £ appear unequivocally related only to the scale 
factor Y c . This result is obtained from Widom's scaling 
law 



{d-l)v = 4> 



(48) 



with d = 3 in our present study. Using Eqs. (45), (47) 
and Widom's scaling law [see Eq. (48)], it is then easy, 
thanks to the formal analogy between the scale dilatation 
method and the analytic hypothesis of the renormaliza- 
tion theory, to validate the well-known universal ampli- 
tude combination previously introduced [see Eq. (38)] 



,(2-1 



= Z a [Z 



d-1 



(49) 



We underline also the microscopic analogy between the 
scale units {/J^ 1 , a c } of the one-component fluid and the 
scale units {ksT c , ai S j ng } used in Monte Carlo simulations 
of the simple cubic Ising-model, where ai s i ng is the spac- 
g lattice size (author?) [73]. Such simulations give 
. — (imiunreoimnrfl v B c , [which compares to Eq. 



o"o = (universal const) x 



((Using)' 



(45)], and £q = (universal const) x aising, [which compare 
to Eq. (47)], to provide a Monte Carlo estimation of the 
above universal ratio (author?) [31, 73]. 

In addition to the universal combination (49), we also 
briefly recall that equivalent universal combinations exist 
between the intcrfacial tension amplitude uq and the heat 
capacity amplitudes A^ as a following form 



1? 



(3c<to (A d 



(50) 



particle energy 



temperature increment 



heat capacity per particle cy,p 

the unique thermodynamic property which can be made 
dimcnsionlcss by only using the "universal" Boltzmann 
factor ks, i.e. without reference to a c and (/3 C ) _1 . There- 
fore, when the singular heat capacity at constant volume, 
normalized per particle, obeys the asymptotic power law 



Ac v ,p 



A ± 

A 0,p 



At* 



1 + 



{lAr-f}] 



(53) 



along the critical isochore, one among the two dimension- 

less amplitudes and is mandatorily a character- 
istic fluid-particle-dependent number. ^From the basic 
hypothesis of the scale dilatation method, this number 
should be related to Y c and Z c in a well-defined manner 
to account for extensive and critical natures of the fluid 
system, as will be shows below [see Eq. (56)]. 

For a 3D Ising-system (author?) [31], the singular 
part of the heat capacity normalized by ks can be ex- 
pressed in units of (ai s i ng ) d . In the case of the one- 
component fluid, the normalized heat capacity is ex- 
pressed in units of (a c ) d , which is the volume of the 
critical interaction cell. For the one-component fluid sub- 
class, the number of particles filling the critical interac- 
tion cell is leading then to define the singular part 
of the master heat capacity per critical interaction cell 
volume (neglecting quantum effects), as follows 



C* s = 



(Y c y z c 



Ac 



V,p 



(54) 



where R+ A « 0.275 ~ (0.537) 3 R~ A (author?) [33, 52, 
53]. Such universal amplitude combinations are related 
to the mixed hyperscaling laws 



a 



d-1 



(51) 



which give common universal features for intcrfacial 
properties (with dimension d — 1) and bulk properties 
(with dimension d = 3). a ~ 0.11 is the universal critical 
exponent associated to the singular heat capacity. The 
above Eq. (51) can be obtained by combining Widom's 
scaling law, 

(d-l)v = <j) 

[see Eq. (48)] and hyperscaling law 

dv = 2 - a (52) 

Equation (52) means that the free energy of a fluctuating 
bulk volume £ d also tends to an universal value [expressed 
in units of (/3 C ) 1 ] for any system belonging to the Ising- 
like universality class. We can then focuss our interest 
in the singular part Acv t p{Ar*) of the heat capacity at 
constant volume expressed per particle, along the crit- 
ical isochore (ignoring the classical background part of 
the total heat capacity at constant volume). Indeed, the 



with Act 



Ac v 



As needed from thermodynam- 



-v, P kn 

ics, such master heat capacity corresponds to a sec- 
ond derivative, C* s = — 9 df-lf ^ of a master free en- 
ergy A* s (T*) with respect to master thermal field T* 
(we do not consider here the critical contribution of 
an additive master constant A* and the regular back- 
ground contribution A* B (T*) characteristic of each one- 
component fluid). Admitting now that the leading 
singular part of the master free energy behaves as 



A 



I T* 



1 + 



{| T *|Aj 



the master 



S ~ q(1-q)(2- 

asymptotic behavior of the heat capacity reads as follows 
(ignoring the critical and background contributions due 
to derivatives) 



C* 8 = 



Z ± 



1 + aZ 



(55) 



The constant values of the master amplitudes are Z^ ~ 

0.1057 ~ 0.537ZI, so that Z A « 0.1967, in conformity 

z + 

with their universal ratio --4- ~ 0.537 (author?) [33] for 

A 

d = 3. Therefore, the corresponding typical asymptotic 
amplitudes A^p in Eq. (53) can be estimated from 



1 A),p 
Z c ks 



(Y c f- a Z± 



(56) 
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where the respective scale factor contributions (i.e., the 
master nature of the critical interaction cell volume char- 
acterized by -j-, and the field scale dilatation along the 
critical isochore characterized by Y c ), are well- identified 
using a single amplitude which de facto characterizes the 
particle. 

Now, for comparison with standard notations used in 
the literature on fluid-related critical phenomena where 
all the thermodynamic potentials are taken per unit vol- 
ume, and not per particle, we also introduce the heat 
capacity at constant volume per unit volume, Acy=i = 
Ac y< p (labeled here with the subscript V = 1). Expressed 

in our above unit length scale [Eq. (17)], the associated 
dimcnsionless form is 



the well-known universal quantities 



a* s x (ty 



(T*) 2 xC* s x (£*f 



zA z iY d 



(61) 



Ac 



Ac 



v=\ 



V,p 



Ac v= 



-d 



k B (a c ) 



-d 



(57) 



where i?+ ~ 0.2696 and w 0.1692, for d = 3 (au- 
thor?) [33]. 

Summarizing the above results for (seven) singular be- 
haviors [surface tension, (±)-correlation length, (±)-heat 
capacity, and (±)-isothermal susceptibility], we note that 
the (five) Eqs. (44), (45), (46), (47), and (58) close the 
hyperscaling universal features along the critical isochore 
above and below T c , in conformity with the two-scale- 
factor universality. Therefore, among the three universal 



exponents v, <j>, and a, only one is readily independent 
Obviously, Ac* v=1 is strictly identical to the usual di- [ scc Eqs _ ( 48 ) and ( 52 )]. The related master - physical 



mensionless form Ac 



heat capacity ACy = NAcy,p of the constant total fluid 
volume V i filled with the constant amount N of parti- 
cles. Using the total Hclmholtz free energy A (T, V, N) 
where T, V, N are the selected three natural variables , 



of the total singular 



we obtain % 



\dTi J 



V,N 



From the corrcspond- 



i.i -i A(T,l,n=&) 

mg quantities normalized per unit volume, — - — y — — . 



we obtain 



£1(4) 

dT 2 



which can be consid- 



ered to define their related singular dimcnsionless parts 



\v)s 



and 



T c Ac 



along the critical isochore. Admitting now that the lead- 
ing singular term of the free energy divergence behaves 



a.s 



(A\* 



1 



0{|Ar*| A }] (ig 



a(l-a)(2-a) I I 

noring critical constant and regular background terms), 
the asymptotic behavior of the singular heat capacity is 



A ± 



1 + aAf | Ar" 



A 



(58) 



Therefore, the leading amplitudes A^ can be estimated 
from 



amplitudes Z a - <jo, Z^ - and Z^ - A^, depend on 
uniquely Y c [see Eqs. (45), (47) and (59)]. 

As a partial but essential conclusion, amplitude cto of 
the interfacial tension, is only characterized by the scale 
factor Y c accounting for the nonuniversal microscopic na- 
ture of each fluid crossing its critical point along the crit- 
ical isochore. 



3.2. The order parameter density dependence 
characterized by the Z c scale factor 

The use of Eq. (35) to estimate the effective parame- 
ters of Eqs. (1) or (10), leads to the Ising-likc expressions 
for the parachor exponent, 



and the (asymptotical) parachor, 



(62) 



Po 
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(63) 



(YcY 



(59) 



which relates A^ only to the single scale factor Y c . How- 
ever the implicit role [see Eq. (56)] of the particle num- 
ber A- filling the critical interaction cell, cannot be ig- 
nored for basic understanding of the master thcrmophys- 
ical properties of the one-component fluid subclass. 

As expected, using Eqs. (45) and (59), we retrieve the 
universal amplitude combinations of Eqs. (50), such as 
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Simultaneously, using Eqs. (47) and (59), we also retrieve 



In spite of the complex combination of scale factors, we 
note that Y c does not appear in the right hand side of 
Eq. (63). However, as clearly shown in the above 3.1, 
Y c is the characteristic scale factor of the critical isochoric 
path where interfacial properties are defined. This amaz- 
ing and important result is entirely due to the hyperscal- 
ing universal feature associated to the critical isothermal 
path, as will be discussed below. 

Indeed, to close the discussion on hyperscaling in crit- 
ical phenomena, we need to introduce the two universal 
exponents i] ps 0.035 and 5 « 4.82, characterizing univer- 
sal features of correlation function at the critical point 
and thermodynamic function along the critical isotherm, 
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respectively (author?) [54]. 77 and 5 are related by the 
hypcrscaling law 



2-77 



1 



1 



(64) 



Equation (64), added to the previous Eq. (52), relate in 
a unequivocal manner the two (independent) exponents 
5 and a describing the thermodynamics, and the two 
(independent) exponents 77 and v describing the correla- 
tions, via d uniquely. In addition, each {thcrmodynamic- 
correlation} exponent pair, either {5; 77}, or {a; v}, char- 
acterizes each independent thermodynamic path to reach 
the critical point, either the critical isothermal line and 
the critical point itself, or the critical isochoric line, re- 
spectively (author?) [27]. 

Obviously, universal values of the corresponding am- 
plitude combinations have been theoretically estimated 
(author?) [31]. We have already given the universal am- 
plitude combination f-R^l = (£ ± ) d (valid along the 
critical isochore), associated to the scaling law dv = 2 — a. 
So that, to close the presentation of the two-scale-factor 
universality, we can also consider the universal amplitude 
combination 



R D = D c p (D p 



(65) 



along the critical isotherm, associated to the hypcrscaling 
law of Eq. (64). Here, we have anticipated (see below) 
the introduction of the leading amplitudes D p (-D£) and 

Dp (D n ) associated to the singular shape of the order- 
ing field along the critical isotherm and to the singu- 
lar decreasing of the correlation function at the critical 
point, respectively. In the D^,, D x notations, i) super- 
script y = c recalls for the non-zero value of the order 
parameter density in a fluid maintained at constant crit- 
ical temperature; ii) decorated hat recalls for the infi- 
nite size of the order parameter density fluctuations in 
a critical fluid maintained exactly at the critical point; 
iii) subscript x = p recalls for a thermodynamic poten- 
tial which is normalized per volume unit and a definition 
of the order parameter density related to the mass den- 
sity, namely A/5 = p ~ Pc (then D p = D, where D is the 
customary notation of this leading amplitude); and iv) 
subscript x = n recalls for a thermodynamic potential 
which is normalized per particle and a definition of the 
order parameter density related to the number density, 
namely Am* = (a c ) d (n — n c ). 

Considering the {thermodynamic-correlation} pairs 

j-Dp;Dpj, or |d^;D„| defined along the criti- 
cal isotherm and at the critical point, and the 
{thermodynamic-correlation} pair {A^;^} defined 
along the critical isochore, from Eqs. (61) and (65), 
a single amplitude characterizes each thermodynamic 
path crossing the critical point, either at constant crit- 
ical temperature, or at constant critical density, respec- 
tively. Considering the {interfacial-bulk} pairs {(To;^.^} 



and {(Tq',^}, the previous section has shown that Y c 
is precisely the single scale factor of the temperature 
field which characterizes the critical isochore. There- 
fore, in conformity with the two-scale-factor universality, 
we are now concerned by the existence of the equiva- 
lent {interfacial-bulk} pairs, which should involve Ising- 
like leading mplitude of the parachor correlations and 
either D c p or D p . Closing their respective Z c -dependence 
demonstrates then that Z c is precisely the single scale 
factor of the ordering field which characterizes the crit- 
ical isotherm. Obviously, the above Eq. (63) where Pq 
appears only Z c -dcpcndcnt is already in agreement with 
such an universal feature. 

Starting with the scaling law, 



/3(6+l) =2 -a 



(66) 



and using the hyperscaling law of Eq. (64), we can recal- 



culate ir„ 



[see Eq. (62]. We obtain 



5 + 1 



1 



(2-77) 



(67) 



with d — 3. The unequivocal link between the Ising- 
like parachor exponent 7r a and either 7/, or (5, is now 
undoubtedly due to mixed hyperscaling along the criti- 
cal isotherm and at the critical point itself, expliciting 
the respective interface (d — 1) and bulk (d) dimensions. 
This Eq. (67) completes the similar mixed hyperscaling 
link = v = [see Eq. (48)] between interfacial 
exponent <f> and either v, or a, along the critical isochore. 

To find the universal amplitude combinations associ- 
ated with Eq. (67), we need to introduce an unambigu- 
ous definition of the parachor correlations from the cor- 
responding Wegncr-likc expansions expressed in terms of 
the order parameter density. The following (master and 
physical) power laws 



(Mt v r>{l + o[(Mt v )$]} (68) 



D° n (Am| v r»{] 

d; (Ap LV r« {l 



O 



(Am^y) 



(69) 



are more appropriate than Eqs. (1), or (10) in the sense 
where Eq. (68) [or (69)] acts as a two-dimensional equa- 
tion of state for the liquid- vapor interface (along the crit- 
ical isochore) . The dimensionless amplitudes and D p 
are called Ising-like parachors to distinguish them from 
dimensional Pq [see Eq. (63)] called parachor. We re- 
call that Eq. (69) refers to the dimensionless interfacial 
tension, a* = Now, the superscript y = a in D\ 

notations, recalls for the thermodynamic definition of the 
interfacial tension of a non-homogeneous fluid where the 
order parameter density spontaneously is non-zero, along 
the critical isochore. As mentioned above, x = n and 
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x = p reflect the two forms Am* and Ap of the order pa- 
rameter density, leading to the r.h.s. forms of Eq. (69). 



The related independence between D° 

di = {z c y° D° 



and is 



(70) 



Using Eqs. (27) and (29) to compare the leading terms 
of master and physical Eqs. (68) and (69), we obtain 



DZ=(Z c )***Z a = (Z C ) 5 ^Z <7 = {Z C )™Z„ 
D° p =(Z c )^Z a = {Z c )^Z a = (Z c )^Z a 

(71) 

As expected, Eqs. (71) provide unequivocal determina- 
tions of D° and Dp 7 from the scale factor Z c (selecting 
either ir a , or S, or 77, as independent exponent). 

We can define in a similar manner the Z c -depcndence 
of D x and D x introduced through Eq. (65). 

First, the amplitude D x are associated to the singular 
behavior of the ordering field in a three-dimensional ho- 
mogeneous fluid in contact with a particle reservoir, fix- 
ing the non-zero value of the order parameter density, and 
thcrmostated at constant (critical) temperature T = T c . 
In that thermodynamic situation, it is established that 
the (master and physical) ordering fields obey the follow- 
ing power laws 



H* = ±Z 



c H \M*f [l + O \\M*ft]} 



(72) 



Selecting then either S, or 77, as an independent exponent, 
Eqs. (76) relate unequivocally each respective physical 
amplitude D£ or D c p . to the scale factor Z c , 

Second, the amplitude D x is associated to the singular 
behavior of the dimensionlcss spatial correlation function 

G* (At* = 0, Ax* = 0, r* = oc (£) ' 



- 2+v at the cx- 



0, Ax* = 0,r* 

act critical point (r is the direct space position, x = n 
or x — p following the order parameter density choice). 
More precisely, introducing the static structure factor 
Xx (T — T c , x — x c , q), where q is the wavenumber in the 
reciprocal space, such as Xx (T — T c , x — x c , 0) takes the 
same dimension as the corresponding isothermal suscep- 
tibility (see below and reference ) , we define the following 
dimensionlcss singular form of the master and physical 
structure factors 



X* (T* = 0, M* = 0, Q* = q*) oc Z G (Q* 



-.77-2 



X* (At* = 0, Am* = 0, q* = qa c ) oc D n (g*)" 2 
, (At* = 0, Ap = 0, q* = qa c ) cx D p (q* f- 2 



\ 1 



(77) 



(78) 



where Zq is a master constant for the one-component 
fluid subclass. Starting from the isothermal susceptibil- 

ities X T,n = (||) T and X t. p = T , associated to 

the order parameter densities Am* and Ap, respectively, 
it is easy to obtain the following relation between the 
amplitudes of the right hand side of Eq. (78) 



A/4 =±£>£|Am 
Ap p 



±m \Apf 1 1 
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(73) 



where Z^ m 252 is a master value for the one-component 
fluid subclass. As in the above case of dimcnsionless in- 
terfacial tension, the r.h.s. forms of Eqs. (73) refer to 
distinct order parameter densities, Am* and Ap, leading 
to the following Z c -dcpcndence between and D° 



D c n = (Z c ) 5 D c p = (Z c )^ D c 



(74) 



In the case of a critical isothermal fluid, a convenient 
rewriting of the leading term in Eq. (73) is (author?) 

[46] 



(av - V P ,c) 
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(75) 



Using Eqs. (27) and (28) and accounting for dual def- 
initions of the ordering field - order parameter density 
with respect to appropriate free energies, the compari- 
son of the leading terms in Eqs. (72) and (73) gives the 
following results 



Dl = (Z c )^+ 1] Z C H = (Z c )^ Z% 
D< = (Z C )^Z C H =(Z c )^Zf I 



(76) 



D„ 
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Similarly, using Eqs. (77), (78), and adding the rela- 
tions between the master isothermal susceptibility X* = 

( aw* ) an< ^ their associated physical dimcnsionless 
forms (neglecting quantum effects), we obtain the fol- 
lowing relations 



(Z c ) d D 7i 

Zr.Dn - 



= z G 

Zg 



(80) 



Each one of Eqs. (80) gives the expected unequivocal 
link between D n , or D p , and Z c . We note that D n , or 
D p , and Z c are true critical numbers, i.e. dimension- 
less quantities defined at the critical point, exactly. One 
among these critical numbers characterizes the selected 
one-component fluid. Therefore, the above link is "ba- 
sic" because it only depends of the hypothesized linear 
relation between master and physical conjugated (order- 
ing field-order parameter density) variables. Eliminat- 
ing Z c between Eqs. (76) and (80) provides the uni- 
versal amplitude combination of Eq. (65), which closes 
the universal features along the critical isotherm and at 
the exact critical point, in conformity with the two-scale- 
factor universality. Finally, Eqs. (56) and (80) are the 
necessary closure equations which unequivocally relates 
the two (independent) leading amplitudes and the two 
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(independent) scale factors characteristics of each one- 
component fluid, selecting rj and a as two (independent) 
critical exponents. 

^From Eq. (71) and Eqs. (76) or (80), associated with 
hyperscaling law of Eq. (67), it is immediate to construct 
the following new combinations between interfacial am- 
plitudes and bulk amplitudes, whose values are expected 
to be universal 
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d-2 + v 
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(81) 



(82) 



To complete our understanding of the universal fea- 
tures related to a (constrained or spontaneous) non-zero 
value of the order parameter density, we must compare 
also the bulk properties of each (liquid-like or gas-like) 
single phase at critical temperature T = T c , and the 
bulk properties of each (liquid or gas) coexisting phase 
in the non-homogenous domain T < T c (admitting then 
a symmetrical one-component fluid close to the criti- 
cal point). In these comparable three-dimensional sit- 
uations where the symmetrized order parameter den- 
sity can take the same non-zero value at two different 
temperatures, the existence of universal proportionality 
(in units of (/? c ) _1 ) is expected for the singular bulk 
free energy E* of a homogeneous phase, either main- 
tained at constant (critical) temperature T = T c [i.e., 



bulk free energy 
V 



cx D° p \&p\ 5+1 oc D c n |Am*| UTi j, or at con 
stant (critical) volume at T below T c [i.e bulk free energy 



1 



Vl,v 



OC 



a(i-a)(2-a) |At*| 2 We account then for the ther- 
modynamic constraints for coexisting phases, expressed 
i i 
i ? 



ordering field as a function of the order parameter den- 
sity along the critical isotherm. Equations (81) and (83) 
imply that the Ising-like parachors D p (or D„, equiva- 

lently), can also be expressed in terms of the ratio B jX 

(b ) s+1 

(or — , cquivalently) , eliminating then D p (or D p , 
cquivalently) . As a matter of fact, despite an explicit 

- 3 

Y c dependence in the amplitudes B oc (Z c ) 2 (Y c ) and 

A ± oc (Y c ) 2 ~ a Z^, their ratio B A ± always takes appro- 
priate forms to ensure the disappearance of the y c -scale 
factor, and only reflect hyperscaling attached to the crit- 
ical isotherm, which is characterized by the Z c -scale fac- 
tor, uniquely. As a consequence, we obtain the universal 
combinations 
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Similarly, we note that the amplitude products T^B S 1 
or [Bn) 8-1 are associated to the "cross" scaling laws 



dv 



7(5 + 1) 



P 
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(86) 



which also reflect hyperscaling attached to the critical 



isotherm. Here above, = and 



(z c y 



are 



the leading amplitudes of the singular behavior of Xt,p 
and Xt n> while 7 — 1.24 is the related critical exponent 
[where I^ are the customary notations along the critical 
isochore, see below Eq. (87]. The (physical) dimension- 
less susceptibilities obey the following power laws 
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b„ J ' ana Ior tne umver - The corresponding (two-term) singular behavior of the 
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sal features above and below the critical temperature 
along the critical isochore, expressed by the universal 
ratio ~ 0.537. As a result, we obtain the follow- 
ing universal amplitude combinations [with B p = B and 
B n = (Z c ) 1 B. where B is the customary notation of 
this leading amplitude, see Eq. (4)] 



master susceptibility X* = (Z c ) Xr n = %cXT,p, is 



X* 
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This amplitude combination is related to the " cross" scal- 
ing laws [see Eqs. (51) and (66)] 



dv 



P P 



= 5+1 



(84) 



0.0248 are the mas- 
ter values of leading amplitudes, with universal ra- 
tio fi = £i ~ 4.79 (author?) [43]. Now, 

the explicit Independence B oc (Zc) 1 (Yc) 13 and 
oc (Z c ) x (Y c ) 7 Z^ desappear in their combination 
Y ± B 5 ~ 1 , due to Eq. (86). This latter product reflects hy- 
perscaling attached to the Z c -scale factor of the critical 
isotherm, uniquely. Introducing the universal combina- 
tion 



Hereabove, Eq. (84) combines exponent ratios ^ or 2-^21 ; 
which caracterizc bulk properties expressed as a func- 
tion of the order parameter density in the nonhomogc- 
neous domain, to the exponent 8 which caracterizes the 



Y + B S ~ 1 D 
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to eliminate D p = D or using Eqs. (81) and (89), we 
obtain the following universal combinations which con- 
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tain the Ising-likc parachors 

= (d° p )^t+b s -i = (z a ) ^ z+ {z M ) & - 1 

(90) 

Summarizing the results for the above (seven) singu- 
lar properties (surface tension, order parameter density, 
ordering field, (±)-heat capacity, and (±)-isothermal sus- 
ceptibility), we note that (three) Eqs. (65), (83), (89), 

and (two) universal ratios ^ = 4=-, ^ = = ^5-, 

close the hyperscaling universal features, at the critical 
point, along the critical isotherm, and in the nonhomo- 
geneous domain, in conformity with the two-scale-factor 
universality. Therefore, among the exponents n, 7r a , 5, 
and exponent ratios ^j^, jg, only one is readily inde- 
pendent [see Eqs. (64), (67), (84), and (86)]. The re- 
lated master - physical amplitudes Z a - D%, Z^ - D c x , 
Zq -D x , and the related master - physical combinations 
( Z ^y+i " WTT> and z x { Z m) S ~ X ~ I^B 5 " 1 , are uniquely 
Z c -dependent (sec Eqs. (71), (76), (81), (85), and (90)]. 

As a partial but essential conclusion, the Ising-likc 
parachor of the interfacial tension, expressed as a func- 
tion of the order-parameter density, is only characterized 
by the scale factor Z c proper to account for nonuniversal 
microscopic nature of each fluid at its critical point or 
crossing them along the critical isotherm. 

4. CONCLUSIONS 

In contrast with all previous studies on the parachor 
correlations, the present estimation of the behavior of in- 
terfacial tension as a function of the density difference 
of the coexisting vapor and liquid phases in the criti- 
cal region, is made without adjustable parameter when 
Qmin _ {(J)^ ; a ^ Z C ,Y C } is known for a selected one- 
component fluid. 

The inter facial-hulk universal features of exponent 
pairs, {4>;a} and {<fr;i>}, or amplitude pairs, {cr^A^}, 
and {oo;^}, indicate that the singularities of the sur- 
face tension, the (thermodynamic) heat capacity, and 
the (correlation) length, expressed as a function of the 
temperature field along the critical isochorc are well- 
characterized by a single characteristic scale factor. Us- 
ing the scale dilatation method, we have shown that this 
(fluid-dependent) scale factor is Y c , precisely. Similarly, 
the interfacial-bulk universal features of exponent pairs, 
{7r a ;<5} and {7T a ; 17}, or amplitude pairs, {D^;D^}, and 

indicate that the singularities of the surface 

tension, the (thermodynamic) ordering field and suscep- 
tibility, and the (correlation) length, expressed as a func- 
tion of the order parameter density along the critical 
isotherm and in the nonhomogencous domain, are well- 
characterized by a single characteristic scale factor. Us- 
ing the scale dilatation method, we have also shown that 



this (fluid-dependent) scale factor is precisely Z c . More- 
over, the desappearence of the isochoric scale factor Y c in 
the estimation of the (Ising-like and effective) parachors, 
is here well-understood in terms of hyperscaling. Y c and 
Z c are two independent characteristic numbers. They are 
fundamental points for future developments of parachor 
correlations. Such results must also be accounted for, 
in equations of the saturated vapor pressure curve, the 
enthalpy of formation of the vapor-liquid interface, and, 
more generally, in ancillary equations where adjustable 
parameters can be estimated using a limited number of 
well-defined fluid-dependent quantities including Y c and 
Z c . 

Since the present approach accounts for complete uni- 
versal features of critical phenomena, thanks to the scale 
dilatation method, in the absence of theoretical predic- 
tion for the surface tension, Fig. 4 may also be useful for 
correlating interfacial properties with master equation 
of the correlation length, incorporating a phenomeno- 
logical contribution of the confluent corrections to the 
asymptotic limit here analyzed. As a special mention, 
the complete classical-to-critical crossover predicted from 
the Field Theory framework can be used with an exact 
knowledge of the density domain where the correlation 
length and the interface thickness reach the order of mag- 
nitude of the microscopic molecular interaction. When 
the two-phase fluid properties change from the critical 
point to the triple point in such a controlled situation, the 
introduction of supplementary parameters, either having 
crossover nature (such as the crossover temperature for 
example), or having empirical origin (such as the acentric 
factor, for example), should then be made to descrimi- 
nate the non-universal character proper to each fluid sys- 
tem revealed from Fig. 4, at large values of M.* LV [or 
from Fig. 3 of Ref. (author?) [29], at large values of 
T*]. However, in all cases, any supplementary parameter 
would be used in conformity to the above master singular 
behavior of the one-component fluid subclass, for which 
the two scale factors arc now specified in terms of thermo- 
dynamic continuity across the critical point (author?) 
[26, 74]. 



Appendix A: PARACHOR CORRELATION IN 
THE NON-HOMOGENEOUS DOMAIN 

We have shown that the power law Z a (A^) 71 " [see 

Eq. (36)], where M% h can be estimated as a modified 

theoretical function of T* < 0, provides an asymptotic 

scaling behavior which agrees with the available experi- 

1 1 f I 

mental results in the range -M LV < -Mp/f, ~ 0.04 (see 
Figs. 4a and b). In this range close to the critical point, 
the Ising-like universal features estimated from the mas- 
sive renormalization scheme are then correctly accounted 
for, in conformity with a (dimcnsionlcss) fluid character- 
ization which only uses the two scale factors Y c and Z c . 

To magnify the relative master behavior at large values 
of M.* LV , i.e. beyond the preasymptotic range M* LV > 
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Figure 5: (Color online) lin-lin scale: a) Residuals (expressed 
in %) for the experimental master parachor from reference 
to the asymptotical master parachor calculated by Eq. (36), 
as a function of the master order parameter density; the full 
black curve corresponds to Eq. (A2) (see text); b) Deviation 
(expressed in %) of the experimental residuals estimated in 
part (a), from reference to the full black curve of Eq. (A2). 
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0.04, we have reported the (%)-residuals 
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as a function of M. 
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A4* h in Fig. 5a, using lin- 
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lin scale. We have also reported the values M.^xt> [ see 
Eq. (41)], M {lf} [see Eq. (42)], and Mac [see Eq. 
(43)], which characterize the finite distances to the crit- 



ical point where I 



th 



40, 



3, and 



-th 



1, respec- 



tively (as previously discussed in 2.4). This figure reveals 
the unambiguous non-universal nature of each fluid at 
large distance from the critical point, i.e., when £^ h < i. 
Moreover, the observed master behavior in the extended 
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ymptotic domain M^}/d ^ AC V < M {lf} ~ 0.16 



i.e., 40 > £t h > 3, needs to introduce a master correction 
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Figure 6: (Color online) lin-lin scale: Residuals (expressed in 
%) for the experimental master parachor of eight HFCs and 
HCFCs (see Refs. (author?) [56, 57, 58, 59]), from reference 
to the master parachor calculated by Eq. (All). 



which can reach —(10 — 20)% at the largest extension 
M {lf} ~ 0.16. 

To account for this correction in a quantitative man- 
ner, we have used the following convenient form, 



R-s* = -100,4 S exp 



-By 



1 



(A2) 



which has two adjustable master parameters (Ay, and 
£?e) to control the following main features: 

i) within the Ising-like prcasymptotic domain, this 
function takes zero-value which preserves the Ising-like 
universal features contained in A4f h , only charac- 
terized by Y c and Z c ; 

ii) beyond the Ising-like preasymptotic domain, this 
function must provide an easy control of supplementary 
adjustable parameters needed to recover each real fluid 
behavior at large distance from the critical point. 

Practically, the master values = 0.575243 and 
By = 0.302034 of the adjustable parameters were thus 
obtained by optimizing the fit on the experimental %- 
residuals of the normal fluids Xe, Kr, N2, and O2, as 
illustrated by the corresponding black full curve in Fig. 
5a. For any other selected fluid, it was then easy to esti- 
mate the related "non-universal" deviation, expressed in 



%, i.e., SRy 



100 



1 



as shown in Fig. 5b. 



We observe that, the larger the value of the renormalized 
order parameter density, the more the fluid behaviors dif- 
ferentiate and the "larger" is the deviation (in amplitude 
and shape) from the master crossover behavior. 

The theoretical estimation of these fluid differences, for 
example in the range of the VLE line which includes the 
temperature value T = 0.7 T c where the acentric factor is 
defined, is out of the present understanding of any fluid 
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Fluid 




Sb,e 




B e 


4>e 


O"0,e 


7T a ,e 


So,e 


PS 


E} 






(mm 2 ) 








(mNm' 1 ) 




(ruNm' 1 ) 






HFC-32 


0.935257 


8.88119 


0.330003 


2.00938 


1.26526 


74.0278 


3.833665 


5.098332 


1,75798 


0, 145884 


HFC-125 


0.9220591 


4.985323 


0.3182363 


1.897228 


1.240295 


52.68437 


3, 896303 


4.342498 


5, 014323 


0, 5406810 


HFC- 134a 


0.8912903 


5.830025 


0.3115459 


1.876448 


1.202836 


50.00357 


3.859691 


4.843207 


1,087740 


0, 1497465 


HFC- 152a 


0.904293 


8.575218 


0.316563 


1.903243 


1.220856 


59.05906 


3.856065 


4.936506 


1,043043 


0, 1492084 


HCFC-123 


0.9785701 


6.054669 


0.2967141 


1.795946 


1.275284 


59.07641 


4.298258 


4.769023 


2,923514 


0, 4197302 


HCFC-124 


0.9073108 


5.270590 


0.289472 


1.768259 


1.196783 


51.72997 


4,116035 


4.912059 


2,015201 


0, 3815439 


HCFC-141b 


0.9246801 


7.477146 


0.3048532 


1.794835 


1.229533 


60.5395 


4, 033156 


5.721538 


1,678934 


0, 1952588 


HCFC-142b 


0.9173884 


7.002359 


0.3062231 


1.903243 


1.223612 


55.82877 


3.996039 


5.114294 


1,939504 


0,2315527 



Table III: Fitting values of the effective exponent-amplitude parameters for power law description by Eqs. (A3), (A5), (A6), 
(A7), and (All) of experimental interfacial properties (see Refs. (author?) [56, 57, 58, 59]) of eight HFCs and HCFCs. 



theory. However, we can use a practical approach de- 
tailed in Ref. (author?) [75], where the main objective 
is to recover consistency with the usual description of the 
surface tension at large distance from the critical point 
by the equation 



o = Co,, 



(A3) 



In the above effective power law, ao,e and <f> e are the 
two adjustable parameters which characterize each one- 
component fluid. Indeed, the noticeable result is a quasi- 
constant value of 4> e which is found in the 1.2 — 1.3 range 
for most of the fluids, but which significantly differs from 
the mean-field value 4>mf = |, as mentioned in our in- 
troduction. In this Appendix, we limit the analysis to 
the S g and a data which are tabulated as a function of 
T data in Refs. (author?) [58, 59], for eight HCFCs 
and HFCs listed in Table I (for the complete study see 
Ref. (author?) [75]). As in the large number of experi- 
mental works related to the determination of the surface 
tension, only the squared capillary length data were ef- 
fectively measured from differential capillary rise method 
(DCRM), as a function of temperature. Then the surface 
tension data have been calculated using Eq. (6), where 
the density difference pl — pv was estimated from pub- 
lished ancillary equations for the liquid density and va- 
por density, along the VLE line. Therefore, introducing 
the critical temperature T c and the critical density p c of 
each selected fluid given in Table, we can re-calculate the 
Ap* LV data at each tabulated value of T c — T, using thus 
the equation 



gpcSg 



(A4) 



where g — 9,80665 ms -2 is the Earth's gravitational ac- 
celeration. 

We have fitted the tabulated a data, using Eq. (A3). 
In addition we have also fitted the tabulated S g and re- 
calculated Ap* LV data, using the following equations 



S ,e I At 



* \<Pe 



(A5) 



Ap 



LV 



B P At* 



(A6) 



Thus, we have checked the consistency of the above re- 
sults, by fitting the tabulated a data as a function of the 
recalculated Ap* LV data, at the same tabulated T c — T 
data, using the following "parachor" equation 



a = E , e {Ap* LV f a 



(A7) 



The adjustable values of cr ,e, 4>e, S , e , <Pe, B e , (3 e , S , e , 
and 7r aje , are reported in Table III and permit to validate 

the interrelations f3 e = 4> e — tp e , 7r a-e = ^1 — , B e = 



gpcSo,, 



and So e = 



(Se)' 



As a conclusive remark, it 



appears that the Eqs. (A3) and (A7) are explicit results 
due to the initial use of "power laws" [see Eqs. (A5) 
and (A6)] to fit the "measured" S g data and Ap* LV data, 
respectively, in the same restricted temperature range. 
The interfacial properties S g , pL~pv, and a, of each fluid 
are then characterized by two amplitude-exponent pairs, 
i.e., four adjustable parameters (in addition to the needed 
critical parameters which characterize the liquid-vapor 
critical point). Moreover, in such a self-consistent result 
of the fitting, the relative error-bar on the surface tension 
data can then be readily approximated by the sum of the 
relative error-bars on the squared capillary length and 
coexisting relative density measurements (thus including 
the relative uncertainty on the p c value, generally of the 
order of 1%). 

In the next step, at each tabulated temperature, we 
have estimated: 

1) the renormalized surface tension 



-'tab 



(a c ) d 1 /3 c cr 



(A8) 



using the tabulated a data, and 
2) the master surface tension 

(Z c )? AplyY" {I- 

Ay, exp —By ' 



(Z c )SAp*, 



In 



[z c y Ap 



LV 



(A9) 
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using the recalculated Ap* LV data. The %-deviations 



at each value M. 



(Z c ) 2 Ap* LV , have been fitted using the following power 
law 



100£* f (M* LV ) 



(A10) 



where and pf are two adjustable parameters, given in 
Table III. As a final result, the parachor correlation now 
reads 



-Bn{-j^-MMl v ])] (All) 



1 — Ay, exp 

(M* LV r 

where the "confluent" correction to the leading term con- 



tains two parts: i) a master contribution which is the 
same for all the pure fluids; ii) a non-universal contri- 
bution which is characterized by the exponent-amplitude 
pair pf, £/ for each pure fluid. Figure 6 gives the corre- 



£♦*. 



tab 
S* 



1 



over the tem- 



sponding residuals = 100 

perature range of the VLE line covered by the capillary 
rise measurements. In this figure, we have also reported 
each error-bar contribution of the experimental accu- 
racy of ±0.2 mNm -1 claimed by the authors (ignoring 
the corresponding contribution of the claimed accuracy 
(±20 mK) on temperature measurements). 

We conclude that the above practical approach based 
on the correct master description of the asymptotic Ising- 
like domain for the one-component fluid subclass, pro- 
vides an easy control of the adjustable parameters needed 
to recover the experimental behavior at large distance 
from the critical temperature. 
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